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SECTION  1 


INTRODUCTION 

The  purpose  of  this  program  is  to  study  both  theoretically  and 
experimentally  the  use  of  phase  conjugate  optics  for  information  processing  and 
image  transmission  applications.  The  objectives  are  to  analyze  and  develop  phase 
conjugation  techniques  for  optical  signal  processing,  to  determine  the  usefulness 
of  phase  conjugation  for  restoring  spatial  information  transmitted  over  multimode 
optical  fibers,  and  to  study  real-time  optical  convolution  and  correlation  of  two 
spatially  modulated  input  waves.  The  contract  was  modified  after  its  inception 
to  include  a  task  covering  the  theoretical  and  experimental  study  of  materials 
for  use  in  performing  phase  conjugation  of  spatially  modulated  optical  waves. 

The  transmission  of  a  “picture"  through  a  multimode  fiber  suffers  from 
several  problems,  two  of  which  are  modal  dispersion  and  mode  mixing.  It  is  well 
known  that  the  process  of  phase  conjugation  -  that  is,  passing  a  spatially 
modulated  wave  through  a  fiber,  conjugating,  and  then  retransmitting  the  beam 
back  through  the  same  fiber  -  can  compensate  for  these  effects  and  produce  a 
restored  image. 

Further,  the  ability  of  one  nonlinear  technique,  degenerate  four-wave  mixing 
(DFWM) ,  to  perform  the  multiplication  of  two  (or  three)  input  signals  in  a 
nonlinear  medium  provides  the  capability  of  generating  the  optical  convolution  or 
correlation  of  spatially  modulated  waves. 

The  use  of  a  nonlinear  optical  process  to  generate  the  optical 
correlation/convolution  is  essentially  a  generalized  Vander  Lugt  Holographic 
filter  technique  in  which  a  nonlinear  medium  operating  as  a  four-wave  mixer 
replaces  a  conventional  holographic  film  recording.  The  key  advantage  of  such  a 
system  is  that  the  "hologram"  may  be  written  and  played  back  almost 
instantaneously  (i.e.,  real-time  signal  processing  can  be  performed). 

In  the  first  year  of  this  30-month  program  we  have  concentrated  our  efforts 
on  the  study  of  photorefractive  materials  for  phase  conjugation  and  on  the 
restoration  of  images  transmitted  through  multimode  fibers.  He  have  developed  a 
detailed  model  of  the  photorefractive  effect  which  will  be  used  as  a  basis  for 
comparing  photorefractive  materials  for  nonlinear  phase  conjugation.  Preliminary 
results  of  four  materials  surveyed  indicate  that  KNbO^  and  BaTi03  are  the 
most  promising  for  the  desired  applications.  Our  own  experience  with  BaTiOj 
supports  the  projected  high  performance  of  this  material. 
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In  our  fiber  measurements,  we  have  successfully  demonstrated  the  use  of 
phase  conjugation  for  the  restoration  of  images  transmitted  through  a  1.75-m 
multimode  fiber.  The  nonlinear  material  was  BaTi(>3,  with  a  demonstrated 
nonlinear  reflectivity  of  ~140%.  This  large  phase  conjugate  return  was  useful  in 
providing  large  signal  levels  at  the  image  plane  of  our  system.  The  preliminary 
results  indicate  that  a  resolution  of  5  lines/mm  was  obtained  (compared  with  a 
resolution  without  the  fiber  of  ~20  lines/mm).  Considerable  improvement  is 
expected  with  more  optimized  imaging  in  and  out  of  the  fiber. 


SECTION  2 


TECHNICAL  PROGRESS 

A.  ANALYSIS  OF  MECHANISMS  AND  MATERIALS  FOR  PHASE  CONJUGATION 

A  large  number  of  nonlinear  mechanisms  have  been  proposed  or  demonstrated 
which  have  promise  for  phase  conjugation  of  spatially-modulated  optical  signals 
via  degenerate  four-wave  mixing  (DFWM).  We  have  performed  a  preliminary  survey 
of  these  mechanisms  along  with  the  performance  of  specific  materials;  the  results 
are  given  in  Table  1.  Several  of  the  tabulated  materials  have  the  desirable 
properties  of  requiring  modest  (CW)  pump  power  levels  and  having  response  times 
in  the  millisecond  or  sub-millisecond  region.  It  is  expected  that  most  of  the 
materials  should  operate  at  spatial  frequencies  approaching  1000  lines/mm  (1  pm 
resolution  in  the  material),  and  thus  be  suitable  for  handling  high  spatial 
bandwidth  signals. 

Among  the  mechanisms  and  materials  listed  in  Table  1,  we  have  chosen  to 
devote  our  initial  attention  to  photorefractive  materials.  The  specific 
advantages  of  these  materials  are  that  they  operate  over  a  relatively  wide 
frequency  range  at  modest  CW  power  levels  and  have  response  times  which  approach 
the  required  value  of  1  msec.  TWo  photorefractive  materials  with  particular  high 
promise  are  BaT103  (demonstrated  high  efficiency),  and  KNb03  (fast  response 
with  moderate  efficiency),  and  these  have  been  included  In  Table  1. 

We  have  undertaken  a  detailed  study  of  the  photorefractive  effect,  whorfe 
goal  is  the  determination  of  optimum  materials  and  configurations  for  phase 
conjugation  of  spatial  modulated  waves  via  DFWM.  As  a  first  step,  a  literature 
survey  was  conducted;  a  list  of  all  references  encountered  and  the  materials 
studied  is  given  In  Appendix  C.  During  the  course  of  our  literature  survey,  we 
collected  pertinent  performance  parameters  for  a  number  of  prominent  materials; 
these  are  tabulated  in  Table  2.  The  intercomparison  of  materials  should  be 
considered  only  qualitative,  because  experimental  conditions  varied  among  the 
cited  references,  and  often  were  not  optimized  for  the  tabulated  parameter.  In 
Table  2,  the  sensitivity  S  is  related  to  the  writing  of  a  holographic  grating, 
and  is  the  change  in  refractive  index  per  unit  of  absorbed  energy.  An 
alternative  parameter  for  characterizing  the  sensitivity  in  a  two-beam 
interference  experiment  is  the  writing  energy  required  to  obtain  1 %  diffraction 
efficiency,  or  the  energy  to  write  a  grating  with  1/e  of  the  saturation 
efficiency.  In  Table  2  wo  have  not  made  a  distinction  between  those  two 
definitions,  although  they  can  lead  to  differences  of  as  much  as  an  order  of 
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magnitude  in  the  tabulated  value.  The  response  time  is  inversely  proportional  to 
the  input  power  (as  discussed  below),  and  is  the  time  required  to  write  a  grating 
with  1/e  of  its  saturated  efficiency.  There  have  been  only  a  few  experiments  on 
DFWM  in  photoref eractive  materials;  the  observed  signal-to-probe  reflectivity  is 
also  tabulated  in  Table  2.  Finally,  we  have  listed  the  relative  importance  of 
the  three  transport  mechanisms  which  participate  in  forming  the  phase  grating  and 
which  determine  the  relative  spatial  phase  with  respect  to  the  interfering 
optical  fields.  These  mechanisms  are  all  accounted  for  in  the  model  described 
below . 

Under  subcontract  to  HRL  on  this  program.  Dr.  Marshal  Sparks  has  developed  a 
detailed  model  of  the  photoref ractive  effect  whose  purpose  is  to  suggest 
experiments  for  the  measurement  of  critical  parameters,  and  to  identify  preferred 
materials  (and  optimum  experimental  conditions).  At  this  stage  the  model 
considers  the  interference  of  two  beams  to  write  a  simple  holographic  grating 
(see  Figure  1).  The  presence  of  a  third  "readout"  beam  simultaneous  with  the 
"writing"  beams  (as  is  present  in  a  typical  DFWM  experiment),  and  the  dynamic 
influence  of  these  beams  on  the  grating  has  not  been  considered.  The  model  in 
its  present  form  can  provide  information  on  the  writ  ng  energy,  response  time  and 
spatial  resolution  for  a  given  system;  further  refinement  is  required  to  model 
the  nonlinear  reflectivity  in  a  typical  DFWM  experiment.  A  preliminary 
manuscript  of  Dr.  Sparks’  written  contribution  is  given  ^.n  Appendices  A  and  B. 

The  theory  which  has  been  developed  provides  both  transient  and  steady  state 
solutions  to  the  system  equations  for  diffusion,  drift  and  the  photo-voltaic 
effect  acting  individually  or  simultaneously.  In  the  model,  the  photoref ractive 
material  is  an  insulator  with  Nq  optical  levels  per  cubic  centimeter  in  the 
electronic  bandgap.  The  levels  are  sufficiently  deep  that  they  are  not  thermally 
excited  into  the  conduction  band.  The  levels  are  called  optical-donor  levels 
because  electrons  from  these  levels  are  photoexcited  into  the  conduction  band. 

The  density  of  the  optical-donor  levels  that  are  filled  with  electrons  is  Ne 
and  the  density  that  are  unfilled  is  N+.  For  example,  Ne  could  be  the  Fe2+ 
density  in  lithium  nlobate,  N+  the  Fe3+.  In  general,  there  are  many  filled 
levels  Ne  and  many  empty  levels  N+.  The  ratio  Ng/N^.  can  be  varied  by 
oxidizing  or  reducing  the  sample  or  by  doping  the  compensating  impurities.  In 
order  to  make  the  crystal  have  charge  neutrality  when  growing,  there  must  be 
negative  impurities  or  imperfections  present  to  compensate  for  the  positive 
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Figure  1.  Typical  geometry  for  hologram  writing  or 

DFWM  in  a  nonlinear  medium.  In  holographic 
language,  Ej  is  the  reference  writing  beam,  Eo 
la  the  image  writing  Beam,  E3  Is  the  readout 
beam,  and  E4  is  the  signal  beam.  In  the 
language  of  DFWM,  Et  and  E3  are  the  pump 
beams,  E2  is  the  probe  beam,  and  E4  is  the 
signal  beam.  The  parameter  A  Is  the  grating 
spacing. 
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charges  of  Che  Iron.  The  density  of  these  levels  that  compensate  for  the  extra 
positive  charge  of  the  Fe is  equal  to  the  value  H+q  of  Kj.  before  the 
Irradiance  is  turned  on.  Subscript  zeros  denote  the  values  of  parameters  at 
t  =  0,  the  time  at  which  the  irradiance  -.a  turned  on. 

Holes  are  neglected  for  simplicity.  Electrons  are  generated  from  the  filled 
optical-donor  levels  to  the  conduction  band  by  the  irradiance  and  recombine  with 
empty  optical-donor  levels.  Trapping  Is  negected  in  this  model.  We  use  the 
common  photoconductivity  terminology  that  electrons  in  traps  can  be  thermally  re¬ 
excited  into  the  conduction  band.  Thus,  the  capture  of  electrons  by  the  empty 
optical-donor  levels  is  called  recombination,  rather  than  trapping.  This 
recombination  into  the  optical-donor  levels  is  not  to  be  confused  with 
recombination  of  holes  and  electrons  across  the  gap. 

The  emphasis  to  date  has  been  on  lithium  niobate  because  there  is  more 
information  available  on  this  material  than  on  other  materials,  some  of  which  may 
be  better  candidates  for  the  intended  applications.  This  information  is  useful 
in  establishing  the  validity  of  the  theory,  which  will  be  done  to  the  extent 
possible  before  applying  the  theory  to  other  materials  such  as  BaTiC>3  and 
KNb03. 

The  appropriate  system  equations  were  identified  and  solved  in  the 
important  limiting  cases  of  short  drift  distance  and  long  drift  distance  (with 
respect  to  the  grating-spacing  distance.  A)  and  of  short  times  and  steady-state. 

the  conditions  for  short-  and  long-diffusion  distances  are  different  for  the 
short-time-  and  steady-state-  limits. 

Simple  physical  explanations  of  the  results  were  developed.  The  grating 
phase  shifts  for  the  various  combinations  of  drift  and  diffusion  and  of  short- 
and  long-transport  distances  were  explained  without  mathematical  analysis. 

Discrepancies  between  the  key  results  of  Amodei1  and  of  Kukhtarev  and 
coworkers2  and  simple  results  obtained  from  general  principles  were  discovered 
and  resolved.  Amodei's  assumption  that  saturation  and  depletion  are  negligible 
is  inconsistent  with  the  steady-state  solution,  which  is  not  obtained  in  the 
limit  of  negligible  saturation  and  depletion.  Kukhtarev's  solution  for  the 
grating  field  Eg  is  not  valid  for  the  case  of  negligible  diffusion  and  open 
circuit  (for  which  Eg  =  0),  but  is  valid  in  the  presence  of  a  small  amount  of 
diffusion. 
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There  are  several  characteristic  times  in  the  photore tractive  effect  that 
come  out  of  the  theory,  the  most  important  of  which  is  the  dielectric  relaxation 
time  which  controls  the  writing  and  erasing  of  the  holographic  grating 

in  all  cases  (diffusion-,  applied  field-,  and  bulk  photovoltaic  effect- 
controlled).  the  dielectric  relaxation  time  is  given  by  (see  Appendix  A, 
Equation  (3.7)): 


e  e 

J  3  -  3  - 

di  4it<j  4irep<n>  * 


(1) 


where  e  is  the  relative  dielectric  constant,  p  is  the  electron  mobility  and  <n> 
is  the  spatial  average  of  the  conduction  band  electron  density.  A  rate  equation 
solution  for  <n>  gives 


<n> 


0_<I>N 
D  eo 

Y  hu> 
r 


(2) 


where  Op  is  the  optical  cross  section,  <I>  is  the  spatially-average  intensity, 
and  Yr  is  the  recombination  rate.  We  note  that  Yr  Is  given  by 


*  vth  S\0  *  <3> 

where  vt^  is  thermal  velocity  of  the  conduction  band  electrons,  and  S  is 

the  capture  cross  section.  Substituting  Equations  (2)  and  (3)  into  Equation  (1) 

gives 


di 


eh(iiv  .  SN 
_ th  +o 

4ire(JO  <I>N 
D  eo 


(*) 
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Several  Important  conclusions  can  be  drawn  from  Equation  (4).  We  first  note  that 

for  a  given  experiment,  the  response  time  Is  Inversely  proportional  to  the 

average  Intensity,  and  the  product  Is  a  constant  (equal  to  the 

writing  energy  of  the  grating)*  One  Important  means  for  diminishing  the  response 

time  Is  to  chemically  treat  the  crystal  ("reduction")  to  increase  the  ratio 

(Ne0/N+0)  of  filled  to  unfilled  levels.  The  ratio 

Neo^N+o»  not  t*ie  individual  values  of  the  densities,  determines  the 

writing  and  erasing  speed  (for  fairly  wide  ranges  of  Neo  and  Nf^). 

In  comparing  candidate  photorefractive  materials,  Equation  (4)  points  out  the 
importance  of  selecting  materials  with  small  values  of  the  quantity  S/iio^.  One 
difficulty  to  overcome  is  the  lack  of  information  on  the  levels  involved,  as  well 
as  their  optical  capture  cross  sections. 

In  the  experimental  phase  of  our  materials  study,  we  have  carefully 
characterized  HaTiC>3  to  determine  the  suitability  of  this  material  for  use  on 
this  program.  To  this  end  we  have  performed  two-wave  mixing  and  four-wave  mixing 
experiments.  The  two-wave  experiments  have  led  to  the  determination  of  the  sense 
of  the  c-axis  direction,  and  have  shown  a  50%  energy  transfer  from  one  beam  to 
another.  In  these  experiments,  the  angle  0  between  the  grating  k  vector  (kg) 
and  the  positive  c-axis  direction  (see  Figure  2)  was  varied  between  0°  and  20°, 
with  the  highest  energy  transfer  obtained  at  the  largest  angle.  On  the  basis  of 
this  performance,  we  are  still  determining  if  two-wave  mixing  has  applicability 
for  image  amplification  on  this  program. 

The  DFWM  configuration  has  been  used  for  most  of  our  image  transmission 
experiments  in  fibers  and  it  is  expected  to  be  used  for  the 
convolution/correlation  experiments.  To  this  end  we  have  measured  the  DFWM 
reflectivity  as  a  function  of  several  parameters  including  angle  between 
pump/probe,  polarization  of  the  beams,  direction  of  c-axis,  intensity  of  the 
beams,  beam  diameter,  etc.  Because  of  the  slow  time  response  of  the  material 
(~0.1  sec  at  the  power  levels  of  our  experiment)  it  was  necessary  to  perform 
these  experiments  on  a  vibrationally  isolated  table. 
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The  Important  results  of  the  experiments  are  that  the  highest  nonlinear 
reflectivity  Is  obtained  with  the  grating  k.  vector  inclined  to  the  positive  c— 
axis  direction,  as  predicted  by  Feinberg  et.  al.3  In  our  experiments  we  were 
only  able  to  obtain  a  20°  angle,  due  to  the  cut  of  the  BaTi(>3  crystal,  and  its 
large  refractive  index.  Published  calculations  give  an  optimum  angle  of  ~50°, 
although  the  reflectivity  at  20°  is  reduced  by  only  a  factor  of  three.  The 
largest  reflectivity  which  we  have  measured  is  133Z,  which  allows  direct  visual 
observation  of  the  reconstructed  images  at  modest  laser  power  levels. 

B.  TRANSMISSION  OF  IMAGES  THROUGH  FIBERS  USING  PHASE  CONJUGATION 

During  this  portion  of  the  program  we  have  been  investigating  the 
transmission  of  images  through  multimode  fiber  waveguides  using  optical  phase 
conjugation  or  wavefront  reversal  techniques.  The  transmission  of  spatial 
information  through  fibers  is  degraded  due  to  mode  mixing  and  modal  dispersion. 

By  making  use  of  the  wavefront  reversal  property  of  nonlinear  optical  techniques 
we  can  compensate  for  these  effects.  In  our  preliminary  experiments  an  image  was 
transmitted  through  a  fiber,  conjugated,  and  made  to  retraverse  the  same  fiber, 
resulting  in  the  reconstruction  of  the  image  at  the  input  end  of  the  fiber.  For 
this  to  be  a  useful  forward  transmission  technique,  operation  with  two  nearly 
identical  fibers  will  be  required.  The  image  will  be  transmitted  through  the 
first  fiber,  conjugated  and  then  transmitted  through  the  second  fiber  which 
reconstructs  the  image  at  its  output . 

The  preliminary  demonstration  was  accomplished  by  utilizing  a  multimode  (83 
pm  diameter,  1.75  m  long)  step-index  fiber.  The  output  of  an  argon  laser  (5145 
A)  was  transmitted  through  a  portion  of  an  Air  Force  resolution  chart  and  imaged 
at  the  input  end  of  the  fiber.  The  fiber  output  was  directed  into  a  crystal  of 
barium  titanate  (BaTiOj)  in  which  DFWM  was  used  to  generate  the  phase  conjugate 
signal.  The  phase  conjugate  signal  then  retraversed  the  same  fiber  and 
reconstructed  the  image  of  the  resolution  chart.  The  results  of  this  experiment 
are  shown  in  Figure  3.  The  first  photograph  shows  the  input  image  of  an  Air 
Force  resolution  chart.  The  second  photograph  shows  the  output  at  the  end  of  the 
fiber.  The  output  is  typical  of  a  highly  multimode  fiber  where  a  large  number 
of  modes  are  excited.  From  the  output,  it  is  impossible  to  discern  that  any 
spatial  Information  has  been  transmitted  through  the  fiber.  The  reconstructed 
image  after  double  passing  the  fiber  is  shown  in  the  last  photograph. 
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Figure  3.  Demonstration  of  image  transmission 

through  a  fiber,  (a)  Image  of  Air  Force 
resolution  chart  at  input  end  of  fiber. 

(b)  Output  after  a  single  pass  through 
the  fiber,  (c)  Reconstruction  of  resolution 
chart  after  output  has  been  phase  conjugated 
and  retraverses  fiber. 


The  geometry  employed  for  the  four-wave  mixing  experiments  in  this  work  is 
shown  in  Figure  4.  It  consists  of  two  counterpropagating  pumps  (Ef,Eb)  and  a 
probe  wave  (Ep)  (which  is  the  beam  with  the  image  information  imposed  on  it) 
incident  on  the  nonlinear  medium.  These  beams  interact  to  generate  the  phase 
conjugate  signal  (Es).  A  photograph  of  the  actual  apparatus  is  shown  in  Figure 
5. 

The  argon-ion  laser  (5145  A)  used  as  the  illumination  source  has  an 
Intracavity  prism  and  etalon  for  single  longitudinal  mode  operating  giving 
coherence  lengths  on  the  order  of  10  m.  In  order  to  have  useful  interference 
between  the  pump  and  probe  the  coherence  length  of  the  laser  must  be  comparable 
to  the  path  length  differences  encountered  in  the  experiment.  For  our 
experiments  the  probe  path  length  is  longer  than  the  path  length  of  the  pump 
beams  by  at  least  the  length  of  the  fiber,  thus  requiring  a  ~10  m  coherence 
length. 

The  laser  output  is  continually  monitored  by  a  scanning  Fabry-Perot  etalon. 
Because  of  crystal  orientation,  p-polarization  of  the  input  beams  is  required  in 
our  experiments.  Therefore,  a  A/2  plate  and  Faraday  rotator  were  used  to  rotate 
the  plane  of  polarization  of  the  argon  laser.  The  Faraday  rotator  was  also  used 
in  conjunction  with  a  polarizer  to  prevent  any  return  beams  from  entering  the 
laser  and  destroying  single  mode  operation.  The  beam  splitter,  BS2,  picks  off 
the  pump  beams.  A  100  cm  f.l.  lens  was  used  to  reduce  the  pump  beam  size  in  the 
crystal.  The  beam  splitter,  BS3,  picked  off  a  secondary  probe  signal  which  was 
used  to  optimize  the  crystal  orientation.  The  beam  splitter,  BS4,  was  used  to 
generate  the  individual  counterpropagating  pumps.  The  pumps  typically  have  a 
power  ratio  of  2:1.  The  throughput  of  the  beam  splitter,  BS2,  is  used  to  obtain 
the  probe  used  in  the  fiber.  This  probe  beam  was  expanded  by  a  beam  expanding 
telescope  to  illuminate  the  Air  Force  resolution  chart.  Following  the  resolution 
chart  is  a  beam  splitter  BS6.  This  element  is  used  to  pick  off  the  return  signal 
which  forms  an  image  on  a  sheet  of  cardboard.  The  image  plane  is  the  exact  same 
distance  from  the  beam  splitter  as  the  original  Air  Force  chart.  An  f.l.,  7.6  cm 
focal  length  relay  lens  was  used  to  capture  the  light  diffracted  by  the  patterns 
of  the  resolution  chart.  This  diffracted  light  would  otherwise  miss  the  entrance 
aperture  at  the  microscopic  objective.  The  microscope  objectives  are  7X,  0.20 
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Figure  4.  Schematic  of  experiment  used  for  image  transmission  through  fibers. 


. 

► 

Figure  5.  Photograph  of  experimental  apparatus  used  for  image 
transmission  through  fibers. 
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N.A.  The  first  microscope  objective  couples  the  light  into  the  fiber.  The  fibers 
are  held  by  a  unique  Hughes  designed  fiber  holder  and  fiber  assembly  which  has 
five  degrees  of  freedom:  x  -  y  -  z  translation  and  azimuthal  and  polar 
rotations. 

In  our  studies  we  have  used  four  types  of  fibers  which  are  representative  of 
the  various  fiber  designs.  The  first  fiber  is  a  "single"  or  lower  mode  number 
fiber  from  Bell  Laboratories.  It  has  a  step  index  with  a  core  diameter  of  8.5  pm 
and  a  numerical  aperture,  N.A.,  of  0.052.  There  are  several  lengths,  the  longest 
being  60  m.  At  the  wavelengths  of  interest,  ~5000  A,  we  can  propagate  three 
discrete  modes,  of  which  two  are  degenerate.  The  number  of  modes  M  is  given  by  M 
■  k2a2(N.A.)2,  where  a  is  the  core  diameter.  The  second  fiber  is  a  multimode 
graded  index  fiber  supplied  by  ITT.  The  core  diameter  is  50  pm,  the  N.A.  is  0.25 
and  the  length  is  6  m.  Theoretically,  approximately  10**  modes  can  propagate  in 
this  fiber  at  the  wavelengths  of  interest.  The  third  fiber  is  a  Corning  85  pm 
step  index  fiber  (N.A.  »»  0.184)  which  is  1.75  m  in  length.  It  also  can  propagate 
~10l*  modes.  The  fourth  fiber  is  a  multimode  step  index  fiber  supplied  by  Quartz 
Products.  The  core  diameter  is  200  pm  and  the  N.A.  is  0.20.  This  fiber  can 
support  on  the  order  of  10s  modes.  We  have  measured  the  properties  of  fiber 
transmission  at  three  different  wavelengths,  4131  A,  5145  A,  and  6471  A  so  that 
we  can  optimize  both  the  wavelength  and  fiber  for  a  particular  experiment.  As 
mentioned  for  the  initial  image  transmission  studies  we  used  the  Corning  85  pm 
f iber . 

A  second  microscope  objective  (see  Figure  4)  quasicolli mates  the  output  beam 
from  the  fiber.  A  5  cm  focal  length  lens  is  then  used  to  reduce  the  probe  beam 
to  a  diameter  smaller  than  the  diameter  of  the  pump  beams  inside  the  crystal. 

The  pumps  typically  make  an  angle  of  75°  with  respect  to  the  c-axis  and  the  probe, 
an  angle  of  68°.  This  implies  that  the  angle  of  the  grating  normal  is  18°  with 
rspect  to  the  c-axis.  The  highest  reflectivity  is  obtained  when  the  grating 
vector  is  inclined  to  the  positive  c-axis  direction,  as  predicted  by  Feinberg  et. 
al.3  Using  this  geometry  we  have  obtained  reflectivities  in  excess  of  133%. 

These  large  returns  have  enabled  us  to  easily  reconstruct  the  image  and  show 
image  amplification.  At  present  we  have  been  able  to  resolve  better  than  5 
lines/mm  (100  pm  spacing),  see  Figure  6.  (The  resolution  of  the  optical  system 
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without  the  fiber  is  20  lines/mm;  and  we  have  seen  in  a  separate  high  resolution 
set-up  greater  than  114  lines/mm  in  BaTi03  by  itself.)  The  imaging  geometry 
used  for  these  preliminary  demonstrations  is  not  optimized,  and  further  work  will 
be  done  with  considerable  improvement  expected.  In  addition,  effects  such  as 
damping  of  the  higher  order  modes  relative  to  the  lower  order  modes  (which  could 
cause  a  scrambling  or  loss  of  information)  and  polarization  scrambling  effects 
will  be  examined.  In  the  latter  case  typical  commercial  fibers  do  not  preserve 
polarization.  Hence  the  linearly  polarized  input  probe  beam  will  emerge  randomly 
polarized.  For  our  DFWM  experiments  the  pumps  are  linearly  polarized  in  the 
plane  of  the  experiment  and  only  this  component  of  polarization  is  conjugated. 

Because  of  the  possible  systems  applications  it  was  of  interest  to  see  if  a 
fiber  could  be  used  as  the  nonlinear  medium  itself.  The  experimental  arrangement 
for  the  DFWM  experiments  using  fused  silica  fibers  as  the  phase  conjugate  medium 
is  shown  in  Figure  7.  The  forward  pump  and  probe  were  brought  in  parallel  to 
each  other  and  were  symmetrically  displaced  from  the  optic  axis  of  the  first 
lens.  This  lens  was  designed  to  couple  the  forward  pump  and  probe  into  the  low 
order  modes  of  the  fiber.  The  backward  pump  was  coupled  into  the  other  end  of 
the  fiber  by  a  microscope  objective.  The  backward  pump  was  also  aligned  to 
launch  only  low  order  modes.  The  fibers  used  were  typically  5  m  in  length  and 
have  been  discussed  earlier. 

In  performing  these  tests  we  encountered  difficulties  in  determining  if,  in 
fact,  a  DFWM  signal  was  observed  because  of  signal-to-noise  problems  fundamental 
to  a  DFWM  experiment  in  fibers.  There  are  two  basic  noise  sources,  one  arising 
from  the  backward  traveling  pump  and  the  other  due  to  Fresnel  refection  of  the 
probe  from  the  input  end  of  the  fiber.  Both  signals  limit  the  DFWM  reflectivity 
that  can  be  detected  and  are  larger  than  the  expected  signal  based  on  a  simple 
calculation  (see  below).  To  understand  these  noise  sources  refer  to  Figure  7. 

The  experimental  arrangement  allowed  us  in  principle  to  spatially  separate  the 
forward  pump,  backward  pump  and  probe.  The  signal  was  separated  from  the  probe 
path  by  the  signal-beam  splitter.  The  problem  was  to  launch  only  the  lowest 
order  modes.  For  best  coupling  it  was  necessary  for  the  forward  pump  and  probe 
to  be  physically  collinear  with  the  optic  axis  of  the  first  lens.  But  for  good 
phase  conjugate  signal  separation  the  probe  should  be  as  far  from  the  optic  axis 
as  physically  possible,  i.e.,  as  close  to  the  edge  of  the  lens  aperture  as 
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Figure  7. 


Schematic  of  degenerate  four 
fiber  as  nonlinear  medium. 


wave  mixing  experiments  using  silica 
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possible.  In  our  preliminary  work  we  attempted  a  compromise  between  these 
conditions.  Unfortunately,  because  of  mode  coupling  In  the  fiber  there  was 
enough  energy  coupled  from  the  backward  pump  into  higher  order  modes  that  the 
energy  in  these  modes  traveled  back  the  same  path  as  the  signal.  Phase  lock 
detection  was  used  (while  chopping  the  probe),  but  the  scattered  backward  pump 
prevented  a  conclusive  measurement.  The  second  source  was  the  probe  retro- 
reflection  off  the  end  of  the  fiber  and  coupling  optics.  We  were  able  to  reduce 
this  retro-reflection  by  incorporating  a  polarizer-quarter  wave  plate 
combination.  The  linearly  polarized  light  from  the  polarizer  passed  through  the 
X^4  plate  at  45°  with  respect  to  the  ordinary  and  extraordinary  axis.  The 
emergent  light  was  then  circularly  polarized  and  impinged  on  the  fiber  end.  The 
reflected  light  was  also  circularly  polarized.  The  light  then  retraversed  the 
A/4  plate  and  emerged  linearly  polarized  but  now  perpendicular  to  the  input 
polarization  and  was  attenuated  by  the  polarizer.  The  signal  polarization  is  not 
preserved  in  the  fiber  and  emerged  with  a  random  polarization  and  the  component 
parallel  to  the  polarizer  was  passed.  Even  with  this  solution  to  the  retro- 
reflection  problem,  we  were  unable  to  detect  a  DFWM  signal,  which  is  not 
surprising  if  one  does  a  simple  calculation  of  the  expected  conjugate  signal 
strength,  as  given  below. 

For  a  nonresonant  material  the  phase  conjugate  reflectivity  is  given  by 

R  =  tan2  |K| %  ,  (5) 

where 

| K J  (cm'1)  -  (32  x  107  J-)  X - tsas£  (w/cm2  )  (6) 

A  (cm)  n 

For  silica  glass  the  linear  and  nonlinear  indices  of  refraction  are: 


and 

n2  *  1.8  x  10~13  esu, 
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giving 


n2n0 

x  -  7T 


4.3  x  10 


-14 


esu  . 


For  fiber  diameters  of  100  ion  and  for  an  Input  power  of  1  W,  1-2  x  102.  Thus 
for  a  wavelength  of  0.5  inn  and  a  fiber  length  of  100  cm,  R  -  10“ 5Z.  In  order  to 
Increase  this  reflectivity  one  could  dope  the  fiber,  thus  increasing  the 
nonlinear  susceptibility. 

As  a  result  of  the  large  phase  conjugate  returns  obtained  in  HaTiOj,  we 
performed  a  separate  experiment  of  see  if  image  amplification  could  be  achieved 
using  a  phase  conjugate  oscillator  configuration.  These  experlements  were  done 
on  an  AFOSR  contract  (F49620-80-C-0041)  and  are  included  here  due  to  their 
relevance  to  this  program.  The  resonator  consisted  of  a  phase  conjugate  mirror 
and  a  conventional  mirror,  shown  in  Figure  8(a).  The  oscillation  occurs  when 
noise  between  the  mirror  is  amplified  by  energy  supplied  from  the  pump  waves. 
Uith  the  onset  of  oscillation  we  have  observed  significant  pump  depletion. 
Typical  operating  parameters  were  ~30  nW  out  of  the  oscillator  with  ~30  mW  in 
each  of  the  pump  beams.  This  resonator  is  unique  because  one  can  readily  place 
two  dimensional  patterns  In  the  cavity  without  destroying  the  oscillation.  To 
illustrate  this,  an  Air  Force  resolution  target  was  placed  within  the  cavity. 
Because  the  light  exactly  retraces  itself  after  one  cavity  pass,  the  resolution 
chart  test  pattern  is  Impressed  on  the  output  and  image  amplification  occurs. 

The  output  from  the  oscillator  in  this  configuration  is  shown  in  Figure  8(b). 
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IMAGE  AMPLIFICATION 
OF  AIR  FORCE 
RESOLUTION  CHART 
AT  OUPUT  OF 
PHASE  CONJUGATE 
RESONATOR 


<b) 


Figure  8.  (a)  Phase  conjugate  oscillator  configuration. 

(b)  Output  at  image  plane  3  from  phase  conjugate  oscillator  with 
Air  Force  Resolution  chart  in  cavity. 
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SECTION  3 


CONCLUSIONS 


In  the  first  year  of  this  program  we  have  analyzed  photorefractlve  materials 
for  phase  conjugation  of  spatially  modulated  beams,  and  we  have  studied  the  use 
of  phase  conjugation  for  restoring  images  transmitted  through  multimode  fibers. 

We  have  developed  a  detailed  model  of  the  photorefractlve  effect  which  will  be 
used  as  a  basis  for  comparing  photorefractlve  materials  for  nonlinear  phase 
conjugation.  The  results  of  the  model  highlight  the  importance  of  characterizing 
the  impurities  involved,  and  determining  material  parameters  such  as  the  electron 
mobility  and  cross  sections  associated  with  the  impurity  sites.  One  element  of 
our  program  of  study  for  the  coming  year  will  be  an  effort  to  measure  or  estimate 
such  required  parameters. 

Our  survey  of  mechanisms  for  DFWM  (Table  1)  has  highlighted  several 
nonlinear  effects  (aside  from  the  photorefractlve  effect)  with  desirable 
properties  for  phase  conjugation  of  spatially  modulated  optical  signals.  We  plan 
to  analyze  and  evaluate  experimentally  several  of  these  techniques  in  the  coming 
year. 

We  have  successfully  demonstrated  the  use  of  phase  conjugation  for 
restoring  Images  transmitted  through  a  1.75  m  multimode  fiber.  Preliminary 
results  indicate  that  a  resolution  of  5  lines/mm  was  obtained,  compared  with  a 
resolution  without  the  fiber  of  ~20  lines/mm.  We  are  devoting  considerable 
attention  to  the  characterization  and  optimization  of  our  imaging  optics.  We 
expect  considerable  improvement  in  resolution  in  the  near  future. 
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ABSTRACT 


A  preliminary  analysis  of  photorefractlvlty  that  vrf.ll  be  used  to  understand 
the  effect,  Interpret  experimental  results,  plus  select  and  improve  materials  and 
device  performance  Is  presented.  Electron  transport  by  diffusion,  drift  in  an 
electric  field,  and  the  photovoltaic  effect  are  considered  both  individually  and 
collectively.  Steady-state -and  short-time  results  are  derived.  Simple  Intuitive 
explanations  of  the  response  time,  sensitivity,  and  spatial  phase  of  the 
photorefractlve  electric  field  are  given.  The  response  time  for  writing  and 
reading  a  grating  Is  shown  to  be  the  dielectric-relaxation  time  when  rewriting  a 
new  grating  Is  neglected.  The  time  response  of  lithium  ntobate  can  be  improved 
by  reduction,  which  increases  the  ratio  of  the  densities  of  filled  (Fe^)  to 
empty  (Fe 3+)  levels. 


APPENDIX  A 


SECTION  1 
INTRODUCTION 


The  purpose  of  this  study  of  photorefractive  materials  is  to  develop 
a  theory  of  photorefractivity,  with  emphasis  on  simple  physical 
explanations  of  the  results,  in  order  to  Interpret  experimental  results, 
suggest  experiments  to  be  performed,  Identify  materials  parameters  that 
should  be  measured,  select  and  rate  types  of  materials  and  specific 
material,  and  seek  Improved  performance.  Photorefractivity,  which  is  the 
change  in  the  refractive  index  of  a  material  with  irradiance,  is  analogous 
to  photoconductivity,  which  is  the  change  in  the  conductivity  of  a 
material  with  irradiance.  A  related  effect,  the  photovoltaic  effect,  is 
the  generation  of  a  voltage  by  irradiance. 

This  interim  report  was  prepared  in  order  to  make  the  mathematical 
details  of  the  treatment  of  photorefractivity  available  at  the  earliest 
possible  time  and  to  solicit  comments  on  the  analysis.  The  details  of  the 
mathematical  treatment  of  a  simple  photorefractivity  model  that  we  believe 
contains  the  essential  fractures  required  to  explain  photorefractivity  in 
LiNb03,  KNbOj,  BaTi03,  and  BSO  (Bi^Si-^o)  are  presented. 

The  results  of  the  present  report  will  be  applied  to  these  materials,  and 
others  of  current  Interest  in  a  subsequent  report.  The  emphasis  to  date 
has  been  on  lithium  nlobate  (LiNb03)  because  there  is  more  information 
available  on  this  material  than  on  other  materials,  some  of  which  may  be 
better  candidates  for  the  Intended  applications.  This  information  is 
useful  in  establishing  the  validity  of  the  theory,  which  will  be  done  to 
the  extent  possible  before  applying  the  theory  as  discussed  above. 

In  a  subsequent  report,  the  speed  and  sensitivity  of  the  various 
materials  will  be  estimated,  tables  of  values  of  material  parameters  will 
be  compiled,  materials  will  be  compared,  and  methods  of  improving  the 
performance  will  be  discussed  along  with  the  effects  of  changing  various 
system  parameters.  An  important  feature  of  the  subsequent  report  will  be 
a  table  establishing  the  operating  regime  for  various  materials  under  diverse 
operating  conditions.  Table  1  lists  the  various  limits  of  laser-pulse 


duration,  electron-transport  distance,  and  electron-transport  mechanism. 

Finally,  the  present  results  will  be  related  to  previous  results  In  the 
literature. 

As  an  example  of  laprovlng  performance,  the  most  important  method  for 
Increasing  the  speed  of  lithium  nlobate  Is  to  reduce  (the  opposite  of  oxidizing) 
the  sample,  which  Increases  the  ratio  of  the  densities  of  "filled"  (Ncq,  for 
Fe2*)  to  "empty”  (N+q,  for  Fe3*)  levels.  This  ratio  N^/N+q,  not 
the  Individual  values  of  Ngo,  determines  the  writing-  and  erasing-speed  (for 
fairly  wide  ranges  of  values  of  and  H+q),  as  seen  In  the  expression 
for  the  dielectric-relaxation  frequency  In  Equation  3.7. 

The  dielectric-relaxation  frequency  determines  the  write-  and  erase-times. 

(Rewriting  a  new  grating  during  erasure  Is  known  to  affect  the  erase  time.) 
There  are  several  characteristic  times  In  the  photorefractlve  effect  including 
the  following:  the  recombination  time  of  electrons  and  eiqpty  optical-donor 
levels,  the  characteristic  time  for  the  irradlance  to  generate  free  carriers, 
the  dielectric  ralaxation  frequency,  the  trapping  time,  and  the  thermal  freeing 
time.  The  time  constant  that  controls  the  writing  and  erasing  times  is  the 
dielectric  relaxation  frequency  for  all  cases,  namely,  diffusion-,  applied- 
field-,  and  photo-voltaic  effect-controlled.  It  Is  Important  to  use  the  correct 
carrier  density,  which  Is  the  carrier  density  generated  by  the  wrlte-or  erase- 
beams,  In  the  dielectric  relaxation  frequency.  It  is  assumed  for  now  that  the 
explanation  of  the  difference  between  write-  and  erase-times  and  between  erase 
times  for  different  c-axis  directions  Is  that  the  writing  of  a  "secondary 
grating"  by  the  erase-  and  dlf fracted-beams  is  correct.  This  effect  will  be 
Investigated  In  a  later  report. 


Table  1.  Operating  Conditions  and  Material  Parameters 


Laser-Pulse 

Short-Time 

Steady-State 

Duration,  t 

e<<Tdi 

t>>Tdi 

Electron-Transport 

Short-Transport 

Lon  g-Transport 

Distance,  A 

Ak  «1 

Ak  »1 

g 

g 

Electron-Transport 

Diffusion 

Applied  Field 

Photovoltaic 

Mechanism 

Effect 

The  dielectric  frequency  given  by  Equation  3.7,  the  result  for  the  space 
charge  electric  field  E  that  gives  rise  to  the  grating,  and  the  relation  for  the 
grating  efficiency  for  a  given  E  determine  the  performance  of  the  photorefractive 
material.  This  grating  electric  field  is  not  to  be -confused  with  the  electric 
field  of  the  optical  beams.  The  response  time,  which  is  the  read-or  write-time 
in  the  absence  of  writing  a  secondary  grating  during  erasure  is  independent  of 
the  grating  wave  vector  kg  and  diffusion  distance  A.  By  contrast,  the  electric 
field  E  depends  on  A  and  kg  in  general. 

The  photorefractive  effect  is  as  follows:  Illuminating  a  material  with 
light  generates  free  carriers,  which  are  assumed  to  be  conduction  electrons  for 
simplicity  in  the  present  treatment.  Holes  can  contribute  to  the  photo¬ 
refractive  effect  in  general,  but  are  often  negligible.  The  conduction 
electrons  are  assumed  to  come  from  such  impurity  levels  in  the  electronic  band 
gap  as  Fe2+  levels  in  LiNb03,  but  could,  in  general,  come  from  direct 
excitation  of  electrons  across  the  gap,  by  one-  or  two-photon  absorption,  for 
example.  These  impurity  levels  will  be  called  optical-donor  levels.  The 
electrons  drift,  if  they  experience  an  electric  field  or  diffuse  to  a  new 
location  in  the  crystal  where  they  recombine  with  empty  optical-donor  levels. 

The  space-charge  electric  field  generated  between  the  electrons  in  their  new 
positions  and  the  enpty,  positively  charged  optical-donor  levels  which  remain  in 
their  original  positions  causes  the  refractive  index  to  change  by  the  electro¬ 
optic  effect. 

In  addition  to  drift  in  an  externally  applied  dc  field  and  diffusion,  the 
photovoltaic  effect  can  contribute  to  the  photoinduced  electric  field  E.  The 
components  of  E  resulting  from  diffusion,  Ed,  drift  in  the  net  electric  field 
(E  +  Eapp,'  where  ®app  *a  an  externally  applied  field  assumed  to  be 
spatially  constant),  Edr,  and  the  photovoltaic  field,  Ep, 

E  -  Ed  +  Edr  +  Ep 

are  discussed  below.  A  typical  value  of  E  is  10**  V/cm.  The  resulting  change  in 
the  refractive  index,  which  is  called  the  grating  index,  is 
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"g  "  2 


(1.1) 


where  nr  is  the  refractive  index  and  r  is  the  appropriate  linear  electro-optic 
coefficient,  having  a  typical  value  of  10“8  cm/V,  which  gives 
ng  =  0.5  (2.3) 3  l<r8  104  -  6  x  10“4 

In  such  optical  processing  as  four-wave  mixing  used  for  phase  conjugation, 
the  irradiance  changes  from  one  position  to  another  in  the  crystal.  For  example, 
the  beat  patterns  of  two  overlapping  beams  causes  a  periodic  irradiance.  To  be 
concrete,  we  assume  a  sinusoidal  irradiance  distribution  for  writing  a  grating 


I  -  <I>  <1  +  M  coskgz) 


where  <  >  denotes  the  spatial  average 


(1.2) 


(1.3) 


and 


kg  ■  Air  nj  sin  0/A,  (1.4) 

with  20  the  full  angle  in  the  material  between  the  direction  of  propagation  of 
the  two  beams,  A  is  the  wavelength  of  the  light  in  vacuum,  and  A/nr  is  the 
wavelength  of  the  light  in  the  material. 

The  electric  field  and  the  resulting  index  of  refraction  ng  have  a  spatial 
frequency  component  at  the  grating  wavelength  kg.  This  index  grating  can 
diffract  an  incoming  light  beam  with  efficiency1  (for  a  thick  grating) 


A-4 


(1.5) 


n 


where  a  is  the  optical  absorption  coefficient,  d  is  the  thickness  of  the  sample, 
20  is  the  full  angle  between  the  two  writing  beams,  nr  is  the  index  of 
refraction  of  the  material,  ng  is  the  change  in  the  index  resulting  from  the 
photorefractive  effect,  and  X  is  the  vacuum  wavelength  of  the  writing  beams. 

There  is  general  agreement  in  the  literature  that  this  description  'f  the 
photorefractive  effect  is  correct.  However,  there  is  disagreement  as  to  such 
details  of  the  model  as  which  electron  transport  mechanisms  —  drift,  diffusion, 
or  asymmetric  photogeneration  of  electrons  (photovoltaic  effect)  —  is  operative 
in  a  given  material  and  why  the  grating  efficiency  increases  or  decreases  with 
increasing  pulse  duration  tp  for  a  given  energy  density  Itp,  rather  than 
remaining  constant  as  predicted  by  early  theories. 

The  model  is  discussed  in  Section  2.  In  order  to  make  the  first  analysis  as 
simple  as  possible,  the  first  case  considered  (see  Section  3)  is  that  of  an 
applied  field  only.  Diffusion,  the  photovoltaic  effect,  depletion,  and 
saturation  are  neglected.  The  case  of  diffusion  only  and  of  photovoltaic  only 
are  consider  In  Sections  4  and  5  respectively  in  order  to  determine  the  effects 
of  these  two  transport  processes  acting  individually.  Depletion  and  saturation 
are  considered  in  Section  6.  The  steady-state  grating  electric  field  with  an 
applied  electric  field,  diffusion,  and  the  photoelectric  field  included  is 
derived  In  Section  7.  The  results  of  Kukhtarev  and  coworkers2  and  of  Amodel2 
are  related  to  the  present  work  and  discrepancies  between  the  three  sets  of 
results  are  resolved  in  Section  8.  Simple  pedantic  models  of  the  photovoltaic 
effect  and  of  dielectric  relaxation  are  given  in  Sections  9  and  10.  A  summary  of 
the  mathematical  results  are  given  in  Section  11.  The  Appendix  contains  a  list 
of  symbols  and  typical  values  of  parameters  for  lithium  niobate. 


Electrostatic  units  (esu)  are  used  in  the  equations.  Numerical  results  are 
given  in  both  electrostatic  units  and  rationalized  tries  units.  Conversion  factors 
are  listed  in  the  Appendix  for  the  convenience  of  the  reader. 

Important  results  are  denoted  by  underscored  equation  numbers  and  used  in 
the  summary  (Section  11). 


SECTION  2 


MODEL  OF  PHOTOREFRACTI VITY 


The  analysis  of  even  the  simplest  model  of  photorefractivity  is  complicated 
when  carried  out  without  approximations  because  the  coupled  material's  equations 
and  Maxwell's  equations  are  nonlinear,  coupled,  partial  differential  equations. 

In  order  to  develop  the  physical  understanding  required  for  materials-iraprovement 
programs  and  operation-improvement  programs,  a  simple  model  is  discussed  first, 
then  the  effects  of  the  assumptions  made  in  the  simple  model  are  discussed. 

The  simplest  model  of  photorefractivity  is  as  follows:  The  material  is  an 
insulator  with  ND  optical  levels  per  cubic  centimeter  in  the  electronic 
bandgap.  See  Figure  1.  The  levels  are  sufficiently  deep  that  they  are  not 
thermally  excited  into  the  conduction  band.  The  levels  are  called  optical-donor 
levels  because  electrons  from  these  levels  are  photoexcited  into  the  conduction 
band.  The  density  of  the  optical-donor  levels  that  are  filled  with  electrons  is 
Ne  and  the  density  that  are  unfilled  is  N+.  For  example,  Ne  could  be  the 
Fe2+  density  in  lithium  niobate,N+  and  Re34-  density,  and  Np  ■  Ne  +  N+  the 
total  iron  density  including  both  Fe2+  and  Fe3+-  In  general,  there  are  many 
filled  levels  Ne  and  many  empty  levels  N+.  The  ratio  Ne/N+  can  be  varied 
by  oxidizing  or  reducing  the  sample  or  by  doping  with  coiqpensatlng  inpurlties . 

In  order  to  make  the  crystal  have  charge  neutrality,  there  must  be  negative 
impurities  or  imperfections  present  to  compensate  for  the  positive  charges  of  the 
iron.  The  density  of  these  levels  that  compensate  for  the  extra  positive  charge 
of  the  Fe3+  is  equal  to  the  value  N+q  of  N+  before  the  irradiance  is  turned 
on.  Subscript  zeros  denote  the  values  of  parameters  at  t  ■  0,  i.e.,  the  time  at 
which  the  irradiance  is  turned  on. 

Holes  are  neglected  for  simplicity.  Electrons  are  generated  from  the  filled 
optical-donor  levels  to  the  conduction  band  by  the  irradiance  and  recombine  with 
eaqpty  optical-donor  levels.  Trapping  is  neglected  in  this  first  model.  We  used 
the  common  photoconductivity  terminology  that  electrons  in  traps  can  be  thermally 
re-excited  into  the  conduction  band.  Thus,  the  capture  of  electrons  by  the  empty 
optical-^donor  levels  is  called  recombination,  rather  than  trapping.  This 
recombination  into  the  optical-donor  levels  is  not  to  be  confused  with 
recombination  of  holes  and  electrons  across  the  gap. 


Figure  1.  Schematic  illustration  of  the  photorefractivity  model 
show  the  various  electronic  levels  in  the  bandgap. 


Saturation  and  depletion  of  the  optical-donor  levels  are  neglected  when 
possible*  However,  they  oust  be  included  in  the  analysis  of  the  steady  state 
because  the  correct  steady  state  is  not  attained  in  general  when  saturation  and 
depletion  are  neglected. 

The  simple  model  is  believed  to  be  sufficient  to  explain  the  performance  of 
iron-doped  lithium  nlobate,  and  it  may  be  sufficient  for  the  other  materials  of 
current  interest.  Extentions  of  the  model  and  the  applicability  to  other 
materials  will  be  considered  in  subsequent  studies  and  reports. 

The  number  of  electrons  per  cubic  centimeter  generated  by  the  light  is  equal 
to  the  number  of  photons  per  cubic  centimeter  absorbed,  which  is  equal  to  the 
absorbed  energy  per  cubic  centimeter,  al,  (where  a  is  the  optical  absorption 
coefficient)  divided  by  the  photon  energy  liw: 

(Bn/atJj  -  (3N+/3E)x  -  al/hu 

It  was  assumed  that  one  free  electron  is  generated  for  every  absorbed  photon. 

This  generation  rate  can  be  written  in  several  convenient  forms  as  follows: 

al/hu  -  a  N  /hw  -  ytN  -  <yt>  N  (1  +  Mcos  k  z)  ,  (2.1) 

u  e  i  e  I  e  g 

where  Equation  1.2  was  used  for  I  and 


“  “  °DNe  »  (2.2) 

Yx  -  OpI/ho)  ,  <Yj>  -  0D  <I>/hu)  (2.3) 

where  Ojj  is  the  optical  absorption  cross  section  of  the  optical  donor  levels. 
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The  change  In  the  electron  concentration  resulting  from  movement  of  the 
electrons  away  from  the  generation  region  (and  from  the  motion  of  the  electrons 
In  general)  Is  described  by  the  charge  continuity  equation: 

3p/3t  -  7  *  J  -  -dJ/d*  (2.4) 

where  the  charge  density  Is 

p  -  e(N+  -  N^j  -  n)  -  e(N+A  -  n)  (2.5) 

where  e  *  |e|,  the  absolute  magnitude  of  the  electron  charge;  N+a  =  N+  - 
N+o ;  and  J  Is  the  electrical  current  (-e  times  the  net  rate  at  which 
electrons  cross  a  unit  area  of  surface).  Substituting  Equation  2.5  Into 
Equation  2.4  gives  the  first  central  equation  ~  Equation  2.8  below.  The 
recombination  rate  Is 

(3n/3t)  .  =  -vsNn 

recomb  th  + 

where  the  recombination  frequency  Is 

\  1  (2'6) 

where  s  is  the  cross  section  for  the  capture  of  an  electron  by  an  empty  optical 
donor  level  (having  density  N+),  and  vt^  is  the  average  thermal  velocity  of 
the  electrons,  which  has  the  value: 
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(2.7) 


v..  2  10  ca/sec 


at  room  temperature.  Combining  these  contributions  gives  the  basic  equations: 


3n/3t  -  3N+/3t  4  e_1dJ/dz  (2.8) 

3N+/3t  -  -3Ne/3t  -  YxNe  -  y^/N^  (2.9) 

J  ■  -epnE  4  eD3n/3z  4  pN^Ic  (2.10) 

3E/3*-4Ae_1p  (2.11) 

p  -  e(N+  -  -  n)  -  e(N+&  -  n)  (2.12) 

E  -  E  4-  E  (2.13) 

app  sc 


where  D  Is  the  diffusion  coefficient,  p  Is  the  photovoltaic  coefficient,  c  Is  the 
unit  vector  along  the  photovoltaic  direction,  and  e  “  |e|.  The  contributions  to 
the  current  In  Equation  2.10  are  as  follows:  -epnE  =  “eVdrlftn  is  the 


usual  electrical-conductivity  current  J  ■  (®;  eDdn/dz  is  the  diffusion  current; 
and  pNeIc  is  the  photovoltaic  current,  as  discussed  in  Section  9. 

These  system  Equations  2.8  through  2.13  are  nonlinear  because  they  contain 
products  N+n  and  nE  (where  E  is  a  function  of  n  and  N+).  Thus,  their 
solution  is  nontrivial,  and  the  complexity  of  the  general  results  would  obscure 
simple  interpretation.  The  following  limiting  cases  will  therefore  be 
considered: 

(a)  kgA  «  1  and  kgA  »  1 

(b)  drift,  diffusion,  and  photovoltaic-effect-nlominated 

(c)  short-time  limit  and  steady  state, 

where  A  is  the  electron-transport  distance,  as  discussed  below. 

For  all  cases,  a  first  approximation  which  will  be  relaxed  below  for  special 
cases,  the  equations  are  linearized  by  neglecting  depletion  and  saturation. 
Depletion  is  the  reduction  in  the  absorption  coefficient  a  resulting  from 
depletion  of  the  electrons  in  the  optical-donor  levels,  and  saturation  is  the 
reduction  in  the  recombination  rate  Yr  resulting  from  filling  of  the 
recombination  levels  (which  are  the  empty  optical-donor  levels  H+.  in  the 
present  model)  with  electrons.  Depletion  is  negligible  when  Ne  8  Neg,  and 
saturation  is  negligible  when  N+  2  N+0*  ^  BDre  specific,  define  the 

changes  K+a  and  Ne ^  as 


(2.14) 


and 


.  U  -  N  . 
eA  e  eO 


(2.15) 


Since  3Ne/3t  ■  -3K*./3t  according  to  Equation  2.9 


N 


+A 


-N 


eA 


(2.16) 


Thus,  the  approximations  of  negligible  depletion  and  saturation  are, 
respectively. 


(2.17) 


lNell  ■  KJ  «  » 


+A1 


eO 


(2.18) 


For  negligible  saturation,  the  optical  absorption  coefficient  a  -  Of)Ne 
Is  constant.  Thus, 


a  -  aDNeQ  ,  for  negligible  depletion.  (2.19) 

For  negligible  saturation,  the  recombination  frequency  vt^8N+  is  constant, 
thus 


y  -  V..8N.-  ,  for  negligible  saturation.  (2-20) 

r  tn  +u 
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SECTION  3 


N 


APPLIED-FIELD  DOMINATED 

In  this  section,  diffusion,  the  photovoltaic  effect,  depletion,  and 
saturation  are  neglected.  The  electron-drift  distance  Is 

Ag  -  UTrE  , 

where  tr  =  l/yr  is  the  recombination  tine  and  pE  is  the  drift  velocity.  The 
two  limits  of  short  drift  distance,  kgA<jr  «  *  an^  long  drift  distance 
kgAg  »  1,  will  be  considered. 

The  steady-state  solution  of  Equation  2.9  with  Ne  s  Neg  and  N+  2 
N+o  is 


n  -  (Yj/Yr)NeQ  -  <n>  (1  +  M  cos  kgz)  (3.1) 

where 


<n>  "  «Yi>/Yr)Ne0 


(3-2) 


typical  value 


t  “  1/v.  sNj.  -  0.2  ns/c,(ppo)  -  20  psec  (3.3) 

r  tn  tv)  T 


« 


« 


Y 


1 
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for  vth  *  10  7  cm/sec  (at  room  temperature),  electron-ion  cross  section  for  a 
changed  recombination  center  s  “  10” 14  cm2,  and  recombination  center  density 
N+0  ■  5  x  10 17  cm-3  for  a  lightly  doped  sample  with  concentration  c+  * 

N+0^Natomic  sites  *  ^  ppm.  For  heavily  doped  samples,  Tr  is  even 
shorter  than  the  value  of  20  psec  in  Equation  3.3.  Thus,  for  pulse  durations  and 
doping  concentrations  of  current  interest,  rr  «  t  is  satisfied  throughout  most 
of  the  pulse  duration.  Furthermore,  n  is  given  by  the  steady-state  value  in 
Equation  3.1,  except  at  very  short  times  which  will  be  neglected  here. 

Physically,  for  negligible  electron  drift  distance,  the  electrons  are  trapped 
where  they  are  generated.  Thus,  setting  J  =  0  in  Equation  2.10  as  a  zeroth-order 
approximation  and  solving  for  n  gives  an  exponential  rise  to  the  value  in 
Equation  3.1  with  rise  time  Tr. 

Having  calculated  n  under  the  zeroth-order  approximation,  J  ■  0,  this  value 
of  n  and  the  nonzero  value  of  J  -  -eunE  in  the  next-order  approximation  cause 
N+  to  continue  to  change  after  n  has  reached  its  steady  state.  For  short 
times(l.e.,  before  the  space  charge  field  Esc  builds  up  to  substantially 
change  the  electric  field  in  Equation  2.9  from  its  original  value  Eq  ” 

Eapp),  Equation  2.10  gives: 


J 


2 

-e  vm 


E 

app 


Substituting  this  J  into  Equation  2.4  and  integrating  over  time,  with  p  *  Pq  * 
0  at  t  ■  0,  gives 

p  -  -ew<n>E  k  M  t  sin  k  z  .  (3.4) 

app  g  g 

The  space-charge  electric  field  resulting  from  the  charge  distribution  in 
Equation  3.4,  obtained  by  integrating  Poissons  equation 
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V  *  s  e3E/3z  ■  4  up 


0-5) 


is 


E 

sc 


’di 


t  E  M  cos  k  z 

app  8 


where  Che  dielectric  relaxation  frequency  is 


(3.6) 


Y „  “  4  ir o/e  -  4  ireu<n>/e  -  4  neya„<I>N  ^/ehujv,  sN .  (3.7) 

ai  L)  eu  t  n  +o  - 

One  of  the  laqiortant  results  of  the  present  study  is  that  the  response  time 
(the  read-  and  write-time  in  the  absence  of  writing  of  a  secondary  grating)  is 
the  dielectric  relaxation  time  **  l/y^.  Consideration  of  the  expression 
for  in  Equation  3.7  can  lead  to  estimates  of  the  response  time  and  to 
suggestions  on  how  to  increase  the  speed,  as  discussed  briefly  in  Section  1. 

This  electric  field  is  in  spatial  phase  (E  ~  cos  kgZ)  with  the  irradiance,  in 
agreement  with  the  results  obtained  below  by  physical  arguments. 

It  is  instructive  to  rederive  this  electric-field  distribution  in 
Equation  3.5  by  physical  argument  in  order  to  gain  intuition.  Once  the  steady 
state  value  of  n  is  attained,  electrons  that  are  generated  at  z  by  the  irradiance 
are  trapped  at  *  +  Ag  on  the  average.  The  trapped-change  density  eN+A 
(see  Equation  2.12)  can  therefore  be  written  as  the  sum  of  two  terms,  a  positive 
term  from  the  N+  at  z  where  an  electron  is  generated  and  a  negative  term  of 
equal  absolute  magnitude  but  displaced  by  distance  Ag,  corresponding  to  the 
recombination  of  the  electron  at  z  +  Ag.  In  other  words,  the 


current  causes  a  buildup  of  the  displaced  trapped  electrons  and  a  corresponding 
build  up  of  the  positive  charges  left  behind  (see  Figure  2).  Integrating  the 
generation  equation  (dH^/dt)^  -  ( 3N+A/3t)gen  “  YlNeO 
for  short  times  gives 


(eN+A)gen  “  e<YI>NeOC(1  +  M  cos  kgz)  *  (3*8) 

The  corresponding  negative,  displaced,  recombination  term  is 

(eN+A)rcmb  =  "e<YI>NeOtfl  +  M  cos  kg(z  "  Adr^  *  (3‘9) 

Adding  these  two  terms  and  using  the  Taylor  expansion  cos  kg(z  ~  Ag)  h 
cos  kgZ  +  Agkg  sin  kgZ  gives 

eN+A  -  -e<YI>  k^AgN^tM  cos  kgz  •  (3.10) 

With  p  3  eN+A  and  <n>  ■  (<Yi>/Yr)Neo»  Equation  3.10  gives 

Equation  3.4.  The  rest  of  the  derivation  of  the  space-charge  field  resulting  in 
Equation  3.6  is  the  same  as  above.  The  charge  resulting  from  n  and  the 
corresponding  N+^  is  easily  shown  to  be  given  by  the  p  in  Equation  3.4  times 
xr/t,  which  is  negligible  for  t  »  xr. 

As  the  time  increases,  the  space  charge  field  Esc  continues  to  increase 
according  to  Equation  3.6  until  it  is  no  longer  negligible.  Then  the 
approximation  E  2  Eapp  is  no  longer  valid,  and  Equation  3.5  is  no  longer 
valid.  The  electric  field  for  intermediate  times  is  more  difficult  to  calculate, 
but  the  steady-state  value  is  again  easy  to  calculate.  In  the  steady-state, 

3n/3t  ■  0  and  3N+/3t  ■  0;  thus.  Equation  2.8  gives  3J/3s  ■  0,  or 
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+ 


(b)  COMPONENTS  OF  p 


-  PEAK  SHIFT 

n 


+ 
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+ 


eN.  A=  (eN+A) 

+A  +A  GENERATION 

+  (®N+* ) 

RECOMBINATION 


(C)  TOTALED 


Figure  2.  Schematic  illustration  of  the  positive  and 

negative  charge  distributions  for  the  case  of 
short  diffusion  distance  (A  k  »1)  showing 
the  total  charge  d istribution8p  that  is  spatially 
phase  shifted  with  respect  to  the  irradiance. 
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J  *  constant  =  Jq 


For  the  present  case  of  drift  only,  setting  D  -  p  *  0  in  Equation  2.10  gives 


J  ”  -epnE 


(3.11) 


From  Equations  3.10,  3.2,  and  J  ■  Jq, 


JQ/eu  <n> 

1  +  M  cos  k  z 


(3.12) 


The  divergence  in  E  for  M  *  1  at  cos  kgZ  =  -1  is  unphysical.  It  arises  from 
the  assuiqption  that  the  transport  (drift)  length  is  zero.  For  nonzero  transport 
length,  M  is  replaced  by  a  number  less  than  unity,  as  suggested  by  Amodei2  and 
shown  by  Kukhtarev  and  coworkers.3  The  field  E  is  not  sinusoidal  except  in  the 
limit  of  M  «  1  (where  the  expansion  (1  +  M  coskgZ)-1  2  1  +  M  cos  kgZ  is 
valid). 

An  important  feature  of  the  field  E  (which  gives  rise  to  the  grating)  is 
obvious  when  E  is  expressed  in  terms  of  Jq,  as  in  Equation  3.12,  although  it 
has  not  been  discussed  previously  to  my  knowledge.  For  open  circuit  conditions, 
the  steady-state  field  E  vanishes  because  Jq  ■  0  (no  current  for  an  open 
circuit).  In  order  to  write  a  steady-state  grating  (i.e.,  have  E  *  0)  in  the 
steady-state  with  drift  alone,  there  oust  be  an  external  circuit  through  which  a 
current  can  flow.  Otherwise,  the  faces  of  the  crystal  are  charged  by  the  nonzero 
current  before  the  steady-state  is  reached,  and  the  current  approaches  zero  as 
the  time  approaches  infinity  and  the  charges  on  the  faces  approach  the  value 
required  to  make  Eg  *  -  Eapp  (i.e.,  E  *  0).  At  intermediate  times,  the 
grating  is  nonzero. 


A  rough  approximation  to  the  full  time  dependence  for  this  case  (of 
kgAdr  «  1,  Eapp  dominated,  negligible  depletion  and  saturation,  and 
J  *  0)  is  obtained  as  follows:  The  charging  of  the  faces  of  the  crystal  gives, 
from  Appendix  A: 


E  =*  E  e 
app 


"V 


(3.13) 


where  the  dielectric  relaxation  frequency  is  given  by  Equation  3.7.  Formally 
replacing  E  *  Eapp  Equation  3.6  by  E  in  Equation  3.13  and  including  the 
spatially  constant  space-charge  field  gives 


“Y..  t 

E  =  e  (E  +  Y..t  E  M  cos  k  z) 
app  di  app  g 


(3.14) 


The  maximum  of  the  cos  kgZ  term  occurs  at 


C  "  Tdi  =  1/Ydi 


(3.15) 


and  has  the  value  (writing  exp  (-1)  to  avoid  confusion  with  the  electron  charge): 


E  =*  exp  (-1)  E  M  cos  k  z 
mx  r  app  g 


(3.16) 
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Since  the  dielectric  relaxation  time  is  the  time  at  which  the  steady-state 
is  approached,  as  an  approximation  to  the  full  time  dependence,  we  multiply  the 
short-time  electric  field  Eshort-t  by  exp  and  add  the 

result  to  f2  =  [l  -  expC-y^^t) ]2  times  the  steady-state  value  Egg 


-Y,.t  -Y..t  , 

E  3  E  „  e  di  +  E  (1  -  e  di  )2 
short-t  ss 


(3.17) 


The  square  of  f  was  used  rather  than  f  to  insure  that  E  2  ^short-t 
at  short  times.  This  approximation  in  Equation  3.17  does  not  correspond  to  an 
exact  solution  of  the  systems  equation.  It  is  simply  an  intuitive  approximation 
to  the  full  time  dependence.  It  does  give  the  correct  short-time  field  and 
steady-state  field,  and  it  does  develop  with  time  constant  ■  1/y<h* 

The  value  of  E  at  intermediate  times  may  not  be  highly  accurate. 

With  Esj,ort_t  given  by  Equation  3.6  and  Egg  given  by  Equation 
3.12,  substituting  in  Equation  3.17  gives 


*  Jl  u 

E  S  y  t  e  EM  cos  k  z 
di  app  g 

+  [-(JQ/ey<n>)/(l  +  M  cos  kgz)] 


(1- 


t 


)2 


(3-l«) 
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Next  consider  the  case  of  kgAdr0  »  1,  where  A^  - 

J,Tr®app*  Since  the  electrons  then  drift  a  distance  long  with  respect  to  the 

grating  distance  kg  the  electron  density  is  approximately  constant  spatially 

for  t  >  x  .  Since  the  free-electron  generation  rate  (3n/3t)  »  y  N  and 

gen  I  e 

recombination  rate  (dn/St)^^  ”  Yrn  are  the  same  as  for  the  case  of  kgA^r 

«  1,  the  value  of  this  constant  n  is  Just  the  average  value  of  n  from 
Equation  3.1 


n  2  <n>  =»  (<YI>/lfr)Ne0  .  (3.19) 

The  recombination  then  occurs  evenly  over  all  z  (rather  than  at  a  short  distance 
from  the  generation  region  as  was  the  case  for  kgAdr  «  1).  Figure  2  (for 
kg^dr  ^  1)  *a  then  replaced  by  Figure  3  for  the  case  of 

kgAdrO  ^  and  the  second  term  in  Equation  3.9  is  replaced  by  zero,  cor~ 
responding  to  spreading  the  charge  eN+Ac  cos  kg(z  -  Adr)  evenly 
throughout  the  crystal.  The  resulting  value  of  eN+A  =  (eN^)^  + 

(eN+A>rcmb  from  Equation  3.8  and  the  modified  Equation  3.9  is 

N^Mt  cos  kgz  .  (3.20) 

With  p  2  *N+£  and  Equation  3.5,  this  gives 

Esc  "  4we<V  Ne0(kaE)tM  8ln  k*z  (3.21) 

o  o 
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which  is  written  as 


E  ”  Y..tE»Msin  k  z  v  (3.22) 

sc  di  A  g  - - - 

where 

E.  5  Y  /uk  (3.23) 

A  r  g  - 

is  the  field  at  which  the  drift  velocity  pE^  is  equal  to  kg-1/Tr. 

The  steady-state  value  of  E  for  arbitrary  Jq  is,  from  Equation  3.11  and  J 

**  Jo* 

E  -  -Jg/eim  .  (3.25) 

The  corresponding  steady-state  value  of  the  drift  distance  is 

A_  -  jjEt  -  Jn/eY  n  .  (3.26) 

Ess  r  u  r 

For  kgAjpgg  »  1,  n  is  approximately  constant  as  discussed  above.  Then 
Equation  3.25  gives  E  =  constant;  thus,  there  is  no  grating  (which  requires  a  z 
dependence  of  E)  in  this  order  of  approximation.  Keeping  higher  order  terms  in  n 
gives  a  nonzero  grating. 

For  kgAjpgg  «  1,  n  is  given  by  Equation  3.1  and  there  is  a 
nonzero  grating,  which  is  given  by  Equation  3.26.  Thus,  Equation  3.25  gives  the 
steady-state  field 


E  -  - 


JQ/eii<n> 


l  +  a  M  coa  k  z 
8 


(3.27) 


where 


a  *  1  for  Jrtk  /ey  «  1 

u  g  r 

a  ■  0  for  Jrtk  /ey  »  1 

o  g  r 

Equations  3.17,  3.22,  and  3.27  give 

E  2  Yj.te  E,M  sin  k  z 
di  A  g 

“Y-.t  , 

-  [(J0/ep<n»/(l  +  a  M  cos  kgz)](l  -  e  ai  ¥  .  (3.28) 


SECTION  4 


DIFFUSION  DOMINATED 

As  in  the  previous  case  of  drift  only,  the  charge  density  for  the  case  of 
kgA,j  «  1  (where  A^  is  the  diffusion  distance  discussed  below)  for 
diffusion  only  can  be  obtained  from  physical  arguments  similar  to  those  used 
above  in  the  discussion  of  Figure  3.1.  For  diffusion,  half  of  the  electrons 
diffuse  to  the  left  and  half  diffuse  to  the  right,  in  contrast  to  the  drift  case 
above  in  which  all  electron  drift  to  the  right  in  Figure  2.  Thus,  Equation  3.9 
is  replaced  by 


(eN. .)  -  -e<Y  >N  rtt[l  +  ~  M  cos  k  (z  -  A.) 

+A  rcmb  I  eO  *■  2  g  '  d 

+  1m  cos  k  (z  +  A  )  ] 


(4.1) 


Using  the  Taylor  expansion 


1  2  2 

cos  k  (z  ±  A  )  3  cos  k  z  ±  k  A,  sin  k  z  -  —  k  A  cos  k  z 
g  d  g  gd  g  2  g  d  g 


and  adding  Equations  4.1  and  3.8  gives 


eN  “  ^  e<YT>  N  nk  ^A  ^Mt  cos  k  z 
+A  2  I  eO  g  d  g 


(4.2) 


With 


Ad2  -  Dtr  (4.3) 

where  D  is  the  diffusion  coefficient.  Using  n  given  by  Equations  3.1  and  3.2, 
Equation  4.2  and  p  2  eN+A  give 

p  2  ^  eDk  2<n>Mt  cos  k  z  .  (4.4) 

2  g  g 

The  corresponding  electric  field  is,  from  Equation  3.5 

E-E  -  t  Y..t  E,M  sin  ,  (4.5) 

sc  2  di  d  g  - 

where 

E  s  I* Jv  (4.6) 

g 


is  the  electric  field  for  which  the  drift  velocity  pE^  is  equal  to  the  grating 
distance  kg'1  divided  by  the  diffusion  time  (Dkg2)-1. 

The  steady-state  value  of  E  for  the  present  case  of  kgA^  «  1  and 
diffusion  only  is  also  easy  to  obtain.  For  diffusion  only,  the  current,  which 
must  be  spatially  constant  in  the  steady  state,  is  zero,  J  ■  Jq  ■  0,  because 
there  is  no  applied  field  and  no  surface  charges.  Setting  J  -  0  and  p  -  0  in 
Equation  2.10  gives 


pnE  “  D3n/3z 


(4.7) 


Since  J  ■  0,  n  is  given  by  Equations  3.1  and  3.2.  With  Equation  4.7  this  gives 


E.Msin  k  z 

v  _  -i - 8 _ 

1  +  M  cos  k  z 
g 


(4.8) 


For  M  «  1,  E  5  Ejj  M  sin  kgZ  is  in  spatial  quadrature  with  the  irradiance. 

When  M  is  not  small,  E  is  not  sinusoidal,  just  as  it  was  not  for  the  case  of 
drift.  By  symmetry,  the  Fourier  component  in  spatial  phase  is  zero. 

The  electron  transport  stops  when  the  drift  in  the  space-charge  field 
increases  to  the  value  at  which  this  drift  cancels  the  transport  from  diffusion 
Thus,  the  characteristic  time  for  the  steady  state  to  be  approached  is  just  the 
characteristic  time  for  charge  buildup,  which  is  the  dielectric-relaxation  time 
■  l/fdi*  wliere  Ydi  Is  defined  in  Equation  3.7.  This  time 
constant  can  be  verified  as  follows:  For  convenience,  set  sin  kgZ  “  1  and  cos 
kgZ  •  0  in  Equations  4.5  and  4.8.  Equating  the  short-time  solution  in  Equation 
4.5  to  the  steady-state  value  in  Equation  4.8  and  solving  for  t  gives  one 
definition  of  the  characteristic  time  (that  would  equal  the  exponential  time 
constant  if  the  E  field  increased  exponentially).  This  gives  t  3  ,  which 

is  approximately  equal  to  the  dielectric-relaxation  time  as  expected. 


From  Equal  loan  1.17,  A. 5,  and  A. ft 


E  5  V  •  1dlt  2  Ed  M  Y 

~Y,.t  2 

+  [E.M  sin  k  */(l  +  M  cos  k  z)  ]  (l  -  e  )  .  (A. 9) 

a  g  g 


Next  consider  the  case  of  kgA^  »  1.  In  this  long-diffusion-distance 
Halt  for  short  tines,  the  electron  density  is  constant  and  is  given  by 
Equation  3.17,  just  as  for  the  case  of  kgA^f  »  1  considered  in  Section  3. 

The  electric  field  at  short  tines  is  therefore  given  by  Equation  3.20. 

E  -  E  ■  Y..t  E  M  sin  k  z  ,  (A. 10) 

sc  di  rg  g  ■* - 


where 


as  in  Equation  3.21. 


Y  /Uk 

r 


g 
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The  steady-state  value  of  E  from  Equation  4.7  with  n  5  const  Is  E  *  0.  A 
nonzero  value  of  E  In  the  steady-state  results  from  the  next  order  correction  to 
n.  An  example  of  the  calculation  of  the  steady-state  value  of  E  is  given  In 
Section  8  In  the  discussion  of  the  results  of  Kukhtarev  and  coworkers. 

From  Equations  3.17  and  4.10 


-T.-t: 

E  s  y  t  e  E  M  sin  k  z  +  E. 

di  rg  g  Dss 


-ir.-t  , 

(i  -  «  dl  )2 


(4.11) 


As  an  example  of  Eqss,  from  Section  8  for  the  case  of  N+q  «  Ne0> 
from  Equation  8.10 


E„__  -  -  (4-rreN lrt/ek_)  M  sin  k_z 


(4.12) 


A—  3] 


SECTION  5 


PHOTOVOLTAIC-EFFECT  DOMINATED 


The  photovoltaic  effect  is  discussed  in  Section  9.  It  is  easy  to  show  that  the 
results  in  Section  3  (for  applied  field  only)  are  valid  for  the  photoelectric  effect 
if  Eapp  is  replaced  by  the  steady-state  photovoltaic  field: 

E  *  A  y  /y  -  pN  <I>/eu<n>  (5.1) 
p  p  r  e 


from  Equation  9.8.  This  prescription  of  replacing  Eapp  by  Ep  can  be 
understood  as  follows:  The  central  feature  of  the  Glass5  explanation  of  the 
photovoltaic  effect  is  that  the  electrons  generated  by  the  irradiance  are  emitted 
preferentially  in  a  direction  given  by  a  unit  vector  c.  The  electrons  are  displaced 
through  an  effective  distance  Ap  in  the  £  direction,  as  discussed  in  Section  9. 

Thus,  Equation  3.9  is  valid  if  is  replaced  by  Ap.  The  prescription  of 
replacing  Eapp  by  Ep  in  the  expression  m  nEapp/Yr  is 

equivalent  to  replacing  A^r  by  Ap,  as  seen  by  making  the  replacement  (5.1)  in 
the  Ajr  expression.  Replacing  Eapp  by  Ep  in  Equations  3.4  and  3.5  and 
integrating  Equation  3.5  for  short  times  gives,  for  Apkg  «1, 

E  =  y..t  E  M  cos  k  z.  (5.2) 

di  p  g  - 

This  result  can  also  be  obtained  from  the  continuity  equation  3p/3t  **  9J/3z,  with 
I  -  -eunE  +  pNeI  3  pNgI  for  short  times.  Notice  that  kgAp  «  1  is  always 
satisfied  for  cases  of  interest  (because  Ap  3  0.1  nra  and  kg  "  1/500  nm). 


As«  the  time  Increases  further,  the  space  charge  builds  up  until  J  -  Jq  - 
constant,  which  is  the  condition  for  the  steady  state  as  discussed  above.  Setting 
J  *  Jo  *n  Equation  2.10  (with  0*0),  solving  for  the  steady  state  field,  and 
using 


pN  <I>/ep<n>  ■  A  y  /p  ■  E 
eo  p  r  p 


gives 


E  “  [E  (1  +  M  cos  k  z)  -  JQ/eM<n>]/(l  +  M  cos  k  z).  (5.3) 

PS  S 


The  first  term  T?p  is  independent  of  z.  Thus,  there  is  no  steady~state 
grating  for  open  circuit  conditions  (Jq  =  0),  just  as  for  an  applied  field. 
However,  there  is  a  nonzero  correction  to  this  result  when  diffusion  and/or  an 
applied  field  are  added.  The  <n>  in  the  numerator  in  Equation  5.3,  which  comes 
from  n  3  <n>(l  +  M  cos  k^z)  is  a  zeroth  order  approximation.  Kukhtarev  and 
coworkers  obtained  the  following  first  order  correction  for  the  case  of  kgA  «  1: 


n 


<n>  [l  +  M  (1  -  O  cos  k  z] 

g 


(5.4) 


A- 34 


For  Eann  and  D  included  and  K„A  »  1,  n  is  approximately  constant. 

app  5 

Then,  replacing  1  +  M  cos  kgZ  by  1  in  the  denominator  in  Equation  5.3  gives, 
for  kgA  »  1 

E  s  E  (1+Mcoskz)-  J/ep<n>  (5.9) 

P  8  J  - 

Equations  3.17,  5.2,  and  5.9  give  for  this  case  of  kgA  »  1, 

-Y..t  -Y..t  2 

E  S  y..te  1  E  M  cos  k  z  +  E  Mfcos  k  *)  (l  -  e  )  .  (5.10) 

'di  p  g  P  ^  8  - 

As  an  interesting  side  point,  for  an  applied  field  the  voltage  at  small 
times  is  Eapp,  an<*  t*ie  v°lta8e  decays  to  zero  as  the  space-charge  field 
builds  up  as  time  approaches  infinity.  The  space-charge  field  builds  up  to 
cancel  the  applied-field,  given  a  net  zero  field.  For  the  photovoltaic  case, 
the  electric  field  is  small  at  small  times  and  approaches  Ep  as  the  space- 
charge  field  builds  up  as  time  approaches  infinfity.  (The  field  Ep  is  the 
space-charge  field  generated  by  the  photovoltaic  effect).  The  space-charge  field 
builds  up  and  cancels  the  photovoltaic  current,  giving  a  nonzero  field. 


SECTION  6 


DEPLETION  AND  SATURATION 

The  conditions  for  no  depletion  of  the  electrons  in  the  Ne  levels  (which 
determine  the  optical  absorption  coefficient  0  =  OjjNe)  and  no  saturation  of 
the  available  recombination  levels  N+  (which  determine  the  effective 
recombination  frequency  vt}jSN+)  are  given  in  Equations  2.17  and  2.18, 
respectively.  The  first  condition  required  for  no  depletion  and  no  saturation  is 
that  the  steady-state  value  (<Yio>/Yr)NeO  the  amplitude  of  n,  from 
Equations  3.1  and  3.2,  must  be  sufficiently  small  for  Equations  2.17  and  2.18  to 
be  valid.  By  conservation  of  electrons,  the  number  of  electrons  N6q 
available  to  be  photoexcited  at  t  =  0  must  always  be  equal  to  the  total  number  of 
electrons 

N  A  =*  <N  >  +  <n>  (6.1) 

eO  e 

where  <  >  is  spatial  average,  as  defined  in  Equation  1.3. 

Taking  the  spatial  averages  of  Equations  2.17  and  2.18  gives 

l<N+A>l  ^  N+0  *  for  ne8li£ible  saturation  (6.2) 

and 

|<NeA>|*|<N+A>|<<  for  negligible  depletion.  (6.3) 

From  Equation  6.1  and  <Ne>  -  Nsq  ■  <Ne  -  Neo>  ■  <N6a>, 

<n>  =  <N  .  >  . 

eA 


(6.4) 


Thus,  Equations  3.2  and  6.4  give 


<".*>  -  «KI>/’rr)»e0 


(6.5) 


which  requires 


<Yj>  «  Yr  ,  for  negligible  depletion  (6.6) 


as  seen  by  comparing  Equations  6.3  and  6.5. 

The  corresponding  condition  for  negligible  saturation  is  obtained  by  writing 
Equation  6.5  as 


<".»>  -  •  (6-7) 

Comparison  of  Equations  6.7  and  6.2,  with  |N+^|  -  |nca|  from  Equation 
2.16,  gives 


<Yj>  «  (Neo^N+0^r  *  for  ne8lfgi61e  saturation.  (6.8) 

The  value  of  <I>  corresponding  to  Equation  6.6  is  <I>  «  I8at,  where 

I  y  hu»/o 

sat  r 

-  2  x  102  s_1  (4.0  x  10_l9J)/5  x  10"18  cm2 

-  1.6  x  10 11  W/co2  . 


(6.9) 


This  very  large  typical  value  of  Isat  2  10 11  W/cra2,  compared  with  typical 
cw  values  of  10  W/cra2  and  rather  large  pulsed  values  of  108  W/cm2  (1  J/cm2  and  10 
msec)  shows  that  saturation  and  depletion  are  not  likely  to  be  problems  in  fairly 
heavily  doped  samples  even  under  short-pulse-duratlon  operation.  It  is  possible 
that  NeQ/N+0  in  Equation  6.8  could  be  small,  thus  lowering  the  saturation 
threshold.  However,  in  most  current  high-speed,  real-time  applications  it  is 
desirable  to  make  N+q  small  and  N6q  large,  so  that  saturation  is 
negligible  (Equation  6.8),  if  depletion  is  negligible  (Equation  6.6). 

In  addition  to  these  restrictions  (Equations  6.6  and  6.8)  on  the  maximum 
irradiance  I,  there  are  restrictions  imposed  on  the  values  of  N’s  In  order  that 
the  values  of  p  2  eN^  at  longer  times  (t  »  Tr)  can  be  obtained.  For 
example,  the  maximum  value  of  E  in  Equation  3.16  for  the  case  of  Jq  ■  0  and 
applied  field  only  is,  to  within  a  numerical  factor  of  order  1,  Eapp  M,  and 
the  corresponding  value  of  In+a]  from  Equation  3.10  with  t  ■  and 

Adr  *  WTrE  5  **TrEappM  is 


N  =  ek  E  M2/4ire 
sat  g  app 


(6.10) 


With  Equations  6.2  and  6.3,  this  gives 


Nsat«  N40> 


N 


eO 


(6.11) 


where  Naat  is  defined  in  Equation  6.10. 

According  to  Equations  6.11  and  6.10,  both  N+q  and  Neo  must  be 
greater  than 


N 


a  6 ( 3__x  105  cm~l) (104 V/cm) (1/300)  (SV/V) 


sat  .  .  ,_~10  t  \l/2 

4x  4.  x  10  (ergcm) 


in16  "3 

10  cm 


(6.12) 
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for  Eapp  =  10  kV/cm,  e  =  (2.3) 2,  sin  0  ■  0.5,  X  *  500  nm,  and  M  2  1.  This 
Is  a  rather  small  value  of  Nflat  —  corresponding  to  ~  10^/5  x  1022  ■  2  x 
10“  7  0.2  part  per  million  of  Fe2*’  and  Pe3*.  Large  dopings  are  commonly  used. 
Thus,  the  usual  doping  are  sufficiently  large  to  avoid  the  long-time  saturation 
and  depletion.  However,  for  small  values  of  M,  the  factor  M2  in  Equation  6.10 
can  make  Ngat  small. 

In  the  general  case,  Nsat  must  satisfy  Equation  6.10,  with 


N  =  ek  E  M/4xe 
sat  g  gmx 


(6.13) 


where  Egmx  Is  the  maximum  value  of  the  z  dependent  electric  field  that 
gives  the  grating  amplitude. 
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SECTION  7 


STEADY-STATE  RESULTS 


The  steady  state  has  been  considered  for  Individual  processes  In  the 
previous  sections*  Since  the  value  of  kgA  for  all  processes  affects  the 
Individual  processes.  It  is  convenient  to  consider  the  steady  state  with  all 
processes  included.  The  key  equations  for  determining  the  steady  state  field  E 
are 

E  ■  (eD3n/3t  +  pNeI  -  J^/epn 

E  *  (eD3n/3z  +  pN  <I>M  cos  k  z  +  eu<n><E>)/evm 
e  g 

E  *  (eD3n/3t  +  pNeI)/eyn  ,  open  circuit 

E  ■  (eD3n/3z  +  pN  <I>M  cos  k  z)/epn  ,  short  circuit 

e  g 

E  ■  (eD3n/3t  +  pN  <I>M  cos  k  z  +  ep<n>E  )/epn,  for  <E> 
e  g  app 

3E/3z  -  -  n) 

»+  -  (YI/Tr)(».N+<)/n) 

and 

<N. .  -  r>  -  0  .  (7.4) 

+A 

Equation  7.1  Is  obtained  from  the  steady-state  condition  J  »  Jq  "  constant 
*  f  (z,t).  The  average  field  <E>  from  Equation  7.1a  Is 

<E>  -  (pN  <I>  -  J_)/ep<n>  .  (7.5) 

e  u 

Solving  Equation  7.5  for  Jq  and  substituting  the  result  into  Equation  7.1a 
gives  Equation  7.1b.  For  an  open  circuit,  setting  Jq  *  0  in  Equation  7.1a 
gives  Equation  7.1c.  For  a  short  circuit  <E>  -  0,  Equation  7.1b  reduces  to 
Equation  7. Id.  For  a  voltage  Vapp  -  *Eapp  applied  with  ohmic 
contacts,  <E>  »  Eappk,  and  Equation  7.1b  reduces  to  Equation  7.1e. 


(7.1a) 
(7.1b) 
(7.1c) 
(7. Id) 
E  (7.1e) 

(7.2) 

(7.3) 


Equation  7.2  Is  Equation  2.11  with  a  p  given  by  Equation  2.12.  Equation  7.3  is 
obtained  from  Equation  2.9  with  3N4./3t  ■  0. 

Equation  7.4  is  the  quasi  charge-neutrality  condition,  which  arises  as 
follows:  In  general,  charge  neutrality  means  that  the  charge  density  p  is  zero. 
Then,  integrating  Poisson's  equation  7  •  eE  »  4irp  gives  E  -  constant  for 
the  current  case  in  which  e=  constant.  (It  is  easy  to  show  that  the  spatial 
dependence  of  e  -  ng  corresponding  to  the  spatial  dependence  of  grating  index 
ng  is  negligible  in  the  cases  considered  here).  Quasi-charge  neutrality  in  the 
context  of  index  gratings  means  that  the  constant  term  <p>  in  p  *  <p>  +  pA  mist 
equal  zero  in  order  to  avoid  the  unphysical  term  E  ~  <p>z  in  the  electric  field. 
Thus,  Equation  2.12  and  <p>  -  0  give  Equation  7.4. 

The  case  of  short  steady-state  electron  transport  distance  is  simple.  [The 
condition  for  this  short  transport  distance  will  be  considered  later. 

Kukhtarev's  condition  C  «  1  and  Cg  <<  1  are  sufficient  if  <n>  <  N+o,  hut 
are  not  necessary  in  general.  Physically  I  would  expect  kgAga  «  1  with 
Agg  -  pEgaTr  +  (Dxr)1/2.]  The  electron  concentration  from 
Equation  3.1  and  3.2  is  n  ■  <n>(l  -  M  cos  kgZ).  Substitution  into 
Equation  7.1a  gives 

E  *  (E  CM  cos  k  z  +  E.  M  sin  k  t  -  Jrt/ey<n>)/(l  +  M  cos  k  z)  (7.6) 
n  p  8  d  g  o  g 

where 

E  =  Ay  /u  ■  pN  <I>/ep<n>  ,  E.  -  Dk  /p  .  (7.7) 

p  p  r  e  d  g 

The  term  Ep 

<n>  (1  +  M  cos  kgz),  ns  discussed  in  the  derivation  of  Equation  5.7.  [The 
consistency  of  E^  in  Equation  7.6  with  Equations  7.2  to  7.4  will  be 
demonstrated  below.  This  consistency  establishes  the  conditions  required  for 
n  *  <n>  (1  +  M  cos  kgZ)  to  be  valid]. 

The  case  of  long  steady-state  electron-transport  distance  requires  some 
care.  Por  example,  formally  setting  n  ■  <n>  in  Equation  7.1  would  give 

E  -  E  (1  +  M  cos  k  z)  ,  for  n  -  <n>  .  (7.8) 

P  8 


CM  cos  kgz/(l  +  M  cos  kgZ)  arises  from  the  deviation  of  n  from 
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This  result  is  incorrect  because  the  deviation  in  n  from  <n>  gives  terms  in 

Equation  7.1  that  are  much  greater  than  Ep  in  general,  as  will  now  be  shown. 
Setting 


n  *  <n>  + 


l6  »  N  -  <N  >  +  N  . 

0  e  e  e5 


in  Equation  7.3  and  using  N+A  =*  N+  -  N+q  gives 


/  N+0  \/<YT><Na>\ 

''V  n5/<n>j(  Yr<n>  )(1  +  M  CO»  V>(‘  +  Ne«/<V>  (7.! 


-  N+0  <n>  -  nfi  . 


with  n6  «  <n>,  which  will  be  verified  below,  and 


N  x  *  r<N  >  M  cos  It  z 
e6  e  g 


(7.10) 


(1  +  M  cos  k  z)(l  +  r  M  cos  k  z)  *  1  +  \  r  M2  +  (1  +  r)  M  cos  k  z 


1  2 

+  —  r  M  cos  2k  z  , 
Z  o  * 


(7.11) 


Equation  7.9  becomes 


-  "  ■<»«-">  +  <««  -  ">5 


(7 .12) 


where 


«N  -  n>  -  N  (g/n  -  1  -  n) 


(7.13a) 
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and 


<Yi><M> 

(N+A  "  n)5  "  N+0  “  Yr<n>  f(1  +  r)Mco8  kgz  +  2  r  M  COS  kg*^  <7*13b> 


with 


g  =  «YI>/Yr)(<Ne>/N+0(l  +  ^  r  M2)  .  (7.14) 


The  factor  r  in  Equation  7.10  is  a  measure  of  depletion.  For  r  «  1,  Ne  s 
-  ^e0*  Since  N^g  ■  -N^+  from  Equation  2.16,  assuming  that 
NeA  2  Ne6  and  ^e^  2  Ne0  in  Equation  7.10  gives 


M  2  n^/n^ 


(7.15) 


Setting  -  n>  equal  to  zero,  to  obtain  quasi  charge  neutrality, 

gives 


n  -  (1/4  +  g)1/2  -  1/2  (7.16a) 

.  1/2 

2  8  for  g  »  1  (7.16b) 

5  8  for  g  «  1  (7.16c) 

These  results  can  be  written  as: 

<n>  s  [<V<Ne>N+0  C1  +  2  r  M2)/Yr]l/2  »  N+o 

for  <Yl>  »  N+QYr/ <^e>  (7.17a) 
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<n>  2  -  \  r  M2)  «  N^,  for  <Yr>  «  n^yv/<a>  (7.17b) 


We  consider  only  the  low-1 rrad lance  limit  g  «  1.  Equation  7.14  and  g  <<  1 


<YI>  <<  N-H)Yr/<Ne>  ^  +  I  r  M2)  » 


(7.18) 


In  which 


<n>  -  (<YI>/Yr)<Nfi>  (l  +  \  r  M2) 


(7.19) 


which  agrees  with  Equation  3.2  In  the  limit  of  negligible  depletion,  <CNe>  = 
Neo,  and  r  «  1.  Equations  7.14  and  7.19  give 


<N+A  " 


irH2 

N  Cl  +  TX  Ji_  cog  k  z  +  m  — - - r  coskz  .  (7.20) 

^  1  +  ^  r  M2  *  ^l+JrM2  8 


Substituting  Equation  7.20  Into  Equation  7.2  and  Integrating  gives 


E  Ci.t-rLJi.  9m  k  z  +  8l"  V 

9  l  j.  -A-  t-m ^  8  i 


1  +f  rM‘ 


1  +  |  r  M 


(7.21) 


where 


E  *  4veN._/ek 
9  +0  g 


(7.22) 
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For  negligible  depletion  the  Inequality  r  «  1  is  valid,  i.e., 

N+O  «  Ne0  from  Equation  7.15.  With  N6a  ■  -N+a  from 
Equation  2.16  and  N+a  ■  m  cos  kgZ  from  Equation  7.20  with 
dg  «  N+Az >  the  term  Ne$/<Ee>  in  Equation  7.9  is  negligible  for 
either  N+q  «  N6q  or  small  modulation.  In  this  case  of  r  «  1, 

Equation  7.21  gives 

E  -  E  M  sin  k  z.  for  g,  r  «  1  .  (7.23) 

<1  g  - 

and  Equation  7.14  and  N+  -  N+q  +  N+a  give 

N  ■  (1  +  M  cos  k  z),  for  g,  r  «  1  ,  (7.24) 

+  +u  g 

where  g  «  1  is  the  small-irradlance  limit  of  Equation  7.17b  and  r  «  1  is  the 
negligible-depletion  limit.  Thus,  for  M  «  1  not  satisfied,  depletion  is  not 
negligible  (N+a/N+o  ■  M  cos  kgZ  is  not  much  less  than  unity). 

The  result  in  Equation  7.8  obtained  formally  (and  incorrectly)  by  using 
n  *  <n>  in  Equation  7.1  is  negligible  with  respect  to  the  correct  result  in 
Equation  7.23,  i.e., 


E  /E  -  k  A  ey  MirepN^.  -  3.7  x  10-5  «  1  .  (7.25) 

p  q  g  p  r  40 


From  Equation  5.5,  we  have  C-p  ■  E^/Eq.  Thus,  Eq  «  E,j  for 
Cf  »  1.  With  Equation  7.25  this  gives 


E  <<  E  <<  E,  for  long  transport  distance 

p  q  d 


(7.26) 


The  following  check  on  the  consistency  of  Equations  7.23  and  7.1  Is  not 
essential  to  the  understanding  of  the  control  results.  The  small  deviations  in  n 
from  the  large  term  <n>  are  retained  and  shown  to  be  small  and  to  account  for  the 
result  in  Equation  7.23.  The  field  E  in  Equation  7.23  comes  from  the  first  term 
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in  Equation  7.1.  The  small  deviation  in  n  from  <n>  must  be  retained  in  Equation 
7.1  in  order  to  obtain  Equation  7.23.  First  consider  the  short-circuit  case  of 
Equation  7.1c.  The  second  term  in  Equation  7.1c  is  negligible  as  shown  in  Equation 
7.26.  Setting  n  =  <n>  (1  +  nc  M  cos  kgZ)  in  Equation  7.1c  and  equating  the 
resulting  E  to  the  E  in  Equation  7.23  gives,  with 
Dkg/u  -  Ed, 

E.n  M  sin  k  z  =  E  M  sin  k  z 
a  c  g  q  g 

Solving  for  nc  gives 


n 

c 


E  /E  .  =  4ireNw.li/eDk 
q  d  +U  g 


=  1/ST  «  1 


(7.27) 


where  ^  is  Debye-le'gth  parameter  of  Kukhtarev  and  coworkers. 

For  the  open-circuit  case  of  Equation  7.1c  defining  n  «  <n>(l  +  ncs), 
where  nC8  =  nc  M  cos  kgZ  +  ns  M  sin  kgZ,  and  using  (1  +  nca)-1  2 
1  -  nca  for  ngc  «  1  and  setting  the  resulting  E  equal  to  the  E  in 
Equation  7.23  gives 


E.n  M  sin  k  z  -  E.n  M  cos  k  z  +  E  M  cos  k  z 
dc  gds  gp  g 

-  E  n  M  sin  k  z  ■  E  M  sin  k  z 

P  s  g  q  g 


-  E  n  M  cos  k  z 
PC  g 


Setting  the  sum  of  the  coefficients  of  the  cos  kz  terras  equal  to  zero  gives 


where  Apkg  «  1  is  always  satisfied  and  Aj^g2  »  1  is  satisfied  in  the  long 
transport  limit.  Setting  the  sum  of  the  coefficients  of  the  sin  kgZ  terms  equal 
to  zero  gives 


r 


where  Eq/E^  «  1  from  Equation  7.27,  n8  «  1  from  Equation  7.28,  and 
VEd  «  1  from  Equation  7.28. 
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SECTION  8 


RESULTS  OF  KUKHTAREV  AND  COWORKERS  AND  OF  AMODEI 

Several  results  of  Kukhtarev  and  coworkers3  and  of  Amodel2  are  discussed 
In  order  to  show  the  relation  between  the  present  results  and  the  results  of 
these  previous  important  papers  and  moreover  to  resolve  Inconsistencies  among  the 
three  sets  of  results. 

Kukhtarev  and  coworkers  consider  only  the  steady  state.  We  set  their 
thermal  generation  factor  8  equal  to  zero  because  the  levels  of  Interest  here  are 
too  deep  for  thermal  generation  of  free  carriers  to  be  Important.  We  also 
neglect  the  small  correction  to  M  for  the  case  of  kgA  «  1  even  though  It 
should  be  simple  to  obtain  this  correction  from  the  present  model.  The  relations 
between  Kukhtarev's  notation  and  ours  are  as  follows: 


f  ♦  Y  N  ♦  N 

'i  D  + 

YR  +  Tr/N+0  CoS  (2kxX  +  +  cos  kgz 

T  .  V  5t  .  Ed/Eq 


N 


N. 


A  +0 
uT/|e|+D 

*  (E  +  E  )/E 
E  app  p  q 


Their  condition 


yt  «  V/N-fC  (8,1) 

Is  valid  only  In  samples  having  many  more  filled  levels  than  eiqpty  levels,  that 
is, 


N 


+0 


«  N 


eO 


(8.2) 
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This  can  be  seen  as  follows*  Setting  3N+/3t  *  0  and  Ne 
Equation  2.9  and  solving  for  N+  gives 


"  Np  +  N+o  in 


YjfiL 


N+  "  Vi  +  V/N+o 


2  (VYr>NDN+0/n 


(8.3a) 

(8.3b) 


where  the  approximate  equality  In  Equation  8.2b  is  valid  when  the  Kukhtarev-  and 
coworkers  unequality  in  equation  8.1  is  valid.  Again  setting  3N+/3t  ■  0  in 
Equation  2.9,  but  retaining  Ne,  and  using  Equations  2.14  to  2.16  gives 


n 


(8.4) 


where 


n 


unsat 


( W%  “  <<YI>/Yr)NeO  (l  +  M  COB  kgZ> 


(8.5) 


Equation  8.3b  can  be  rewritten  as 


"unsat (ND/NeO)  1  +  N  /N 


(8.6) 


The  exact  result  in  Equation  8.4  reduces  to  the  approximate  result  in  Equation 
8.6  if  Nd  2  Neg.  Then  N+g  <<  Neg,  the  factor 

N+^/Neo  Is  negligible  in  the  negligible-saturation  limit,  and  the 
factor  Ng/Neg  =  l  in  Equation  8.5. 

Equation  9  of  Kukhtarev  and  coworkers  with  m  ■  M  is  in  agreement  with 
Equations  3.12,  4.8,  and  5.3.  Thus,  our  results  and  theirs  agree  for  the  short- 
transport  distance  limit  of  kgA  «  1.  Their  Equation  9  Is  valid  even  when 
their  approximation  N+g  <<  Ng  is  not  satisfied. 


*•"  c°n8lder  «>•  l‘-lt  V  »  1.  Fro.  Equation  8.6  with  the  Knkhtare. 
and  coworkers  approximation  HD/Ne0  -  1,  i„  order  to  have  n  5  <n>  we  test 
have 


N+A  ’  N+0  M  cos  kgz  (8.7) 

so  that  the  factor  (1  +  N+a/N-K))"1  in  Equation  8.6  cancels  the  factor 
(1  +  M  cos  kgz)  in  n^ga,.  (see  Equation  8.5).  Notice  that  in  this  case 
of  N+q  «  nd,  the  steady  state  is  not  reached  until  N+0  is  "fully 
saturated"  for  M  2  1,  that  is, 

N+  =  N+o  C1  +  M  cos  kgz)  .  (8.8) 

The  charge  density  from  Equation  2.5  is 

P  -  e(N+A  -  n).  (8.9) 

With  n  «  N+a,  substituting  Equations  8.9  and  8.8  into  Equation  2.11  and 
Integrating  gives 


E  =  <E>  -  E  M  sin  k  z 
9  8 


(8.10) 


where  <E>  is  the  obvious  constant  of  integration  and 


E  =  4ireN  /ek 
q  +0  g 


(8.11) 


Is  In  agreement  with  Equation  7  of  Kukhtarev  and  coworkers. 
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Finally,  notice  that  Equation  8.10,  which  is  Equation  7  of  Kukhtarev  and 
coworkers,  is  incorrect  for  the  case  of  an  applied  field  only  and  an  open  circuit 
because  Equation  2.10  with  J  -  D  =  p  =*  0  gives  E  ■  0  for  n  *  0.  The  diffusion 
term  in  Equation  2.10  and  the  small  deviation  ng  in  n  from  <n>  are  needed  to 
establish  the  field  in  Equation  8.10  even  though  D  does  not  appear  in  the  result 
(because  the  1/D  in  ng  cancels  the  explicit  factor  D  in  Equation  8.13). 


Amodei2  obtained  the  diffusion-dominated  result 


E  ,  M  sin  k  z 

E  =  -d _ g _ 

1  +  M  cos  k  z, 
8 


(8.12) 


even  in  the  limit  of  A^^^kg  »  1.  This  result  disagrees  with  the  present 
result  in  Equations  7.21  or  7.23  and  with  the  result  in  Equation  8.10  from 
Kukhtarev  and  coworkers.  The  reason  for  the  disagreement  is  that  Amodei  assumes 
that  Neo  (his  N)  and  N+g  (in  his  t  -  l/yr)  are  constant.  Then 


n 


<n>  (1  +  M  cos  k  z) 
8 


(8.13) 


rather  than  n  *  <n>,  as  obtained  here  by  Kukhtarev  and  coworkers 

This  result  in  Equation  8.12  is  not  correct  in  general  because  the  changes 
in  the  values  of  N+  and  Ne  cannot  be  neglected  in  the  steady  state  in 
general.  For  example,  N+g  was  fully  saturated  in  the  case  above,  according 
to  Equation  8.8.  Equation  8.13  is  physically  implausible  for  this  case  of  kgA 
»  1  because  the  transport  distance  is  greater  than  the  distance  kg  over  which 
<I>  M  cos  kgZ  changes  substantially  and  the  electrons  cannot  remain  in  the 
region  in  which  they  are  generated  if  they  diffuse  or  drift  such  a  long  distance 
before  recombining. 
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The  result  In  Equation  8.12  Is  also  Inconsistent  with  the  value  of  E 
obtained  from  Poissons  equations.  For  N+q  =  constant,  N+q  =  0,  and 
Equation  2.5  gives 


P 


-en  =  -e<n>(l  +  M  cos 


k  z) 
8 


(8.14) 


Substituting  the  z-dependence  term  into  Equation  2.11  and  integrating  gives 


» 


■ 


E  *  -(4ire<n>M/ek  )sin  k  z  , 

8  8 

which  does  not  agree  with  Equation  8.12.  Finally,  the  constant  term  in  p  in 
Equation  8.14  does  not  vanish,  as  required  by  quasi-damage  neutrality. 
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SECTION  9 


PHOTOVOLTAIC  EFFECT 

In  photovoltaic  materials,  a  space-charge  field  is  believed  to  be  generated 
by  asymmetric  photoemission  of  conduction  electrons.  1  In  such  ferroelectric 
materials  as  LiNbOj  and  BaTiOj,  ions  can  have  potentials  that  are  asymmetric 
along  some  direction.  For  example,  the  titanium  ion  in  BaTiOj  is  displaced 
from  the  center  of  the  lattice  cell  in  the  ferro-electric  state.  Thus,  the 
photogenerated  electrons  are  emitted  preferentially  in  one  direction,  which  gives 
rise  to  the  photovoltaic  field  Ep  as  follows:  By  definition  of  the  momentum 
transport  time  x^,  the  Initial  momentum  is  lost  in  time  x^  on  the  average. 

The  displacement  of  the  photo gene rated  electrons  during  the  time  x^  has  a 
typical  value  £p  =  vx^  =*  (107  cm/s)(10-lt*  s)  =*  1  nm  (=  10  A).  All  of  the 
electrons  are  not  photoemitted  along  the  same  direction,  but  there  is  a  greater 
probability  of  emission  along  the  photovoltaic  direction.  Thus,  an  effective 
distance  Ap  is  defined  as 

A  »  p  A  -  p  A  ~  (p  -  p  )  A  =  p,  A  (9.1) 

p  +  p+  -  p-  +  -  p  4—  p 


where  p+  and  p_  are  the  probabilities  of  emission  in  the  +  and  -  directions, 
respectively. 

A  typical  value  of 


A  =  0.08  nm  (=  0.8  A)  (9.2) 

P 


is  consistent  with  typical  observed  photovoltages,  as  shown  below.  Equations  3.1 
and  9.2  give  P+_  »  Ap/ Ap  *  0*08  nra/1  nm  ■  0.08  as  the  relative 
probability  for  emission  in  the  positive  photovoltaic  direction.  The  short 
distance  of  0.08  nm,  which  is  less  than  an  atomic  spacing,  is  not  the  actual 
transport  distance,  which  would  be  unrealistically  short.  Rather,  it  is  the 
weighted  distance,  as  in  Equation  9.1. 


A-55 


In  the  literature,3  It  Is  common  practice  to  model  the  photovoltaic  effect 
by  including  the  term 


J 

~P 


(9.3) 


in  the  current,  where  p^  is  a  constant  and  &  is  a  unit  vector  in  the  opposite 
direction  from  the  preferred  electron  emission  direction.  The  following 
derivation  shows  that  Equation  9.3  is  consistent  with  Glass'  preferential 
emission  model  if  the  constant  p^  is  replaced  by  p  Nde ,  where  p  is  a  constant 
and  N^e  is  the  trapped  electron  density: 


J  »  p  N  ii  (9.4) 

P  e 

Equation  9.4,  but  not  Equation  9.3,  is  equivalent  to  Glass'  expression  J  *  *cal 
because  the  absorption  coefficient  a  is  proportion  to  Ne.  Equations  9.3  and 
9.4  are  equivalent  in  the  limit  of  no  depletion  (Ne  »  NeQ  *  constant). 

Equation  9.4  is  obtained  by  considering  the  generation  rate  of  electron  by 
the  field 


(3n/3t)T  -  (a _/hoi)I(z  -  A  )N  (as  -  A  )  .  (9.5) 

ID  pep 


The  electrons  generated  from  optical-donor  levels  at  z  -  Ap  appear  at  the 
displaced  position  z.  To  Taylor  expand  in  Ne  gives 


(3n/3t)_  -  o  I(z)N  (z)hoj  -  3(A  o  I(z)U  (z)/hui)/3z  , 

l  u  e  p  u  e 


which  is  of  the  usual  form 


3n/3t  *  a  IN  /hu>  -  e  ^3J/3z 
D  e 
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if  we  set 


J  -  J  -  pN  Ic  ,  (9.6) 

~p  e 

where  £  *  -z  and 

p  ■  eA  o  f ho»  .  (9.7 

p  D 

We  shall  use  Jp  in  Equation  9.6,  rather  than  the  previous  expression  in 
Equation  9.3. 

The  photovoltaic  field  that  develops  in  the  steady  state  under  open-circuit 
conditions  (that  is,  J  ■  <n>epE  +  pNg<I>  ”0)  is: 

E  -  pN  <I>/<n>ep  “  A  y  f\i  .  (9.8) 

p  e  p  r 

The  time  constant  for  Ep  to  approach  the  steady  state  value  in  Equation  9.8  is 

the  dielectric  relaxation  frequency  Y<ji,  defined  in  Equation  3.7,  as 

usual. 

A  small  photovoltage 


E  =  4wE<n>A  ft 
psm  p 


<WS 


(9.9) 


is  developed  in  the  short-time  xr,  as  seen  in  the  following  model.  Figure  4 
schematically  illustrates  the  charge  distributions  resulting  from  the 
displacement  of  the  electrons,  denoted  by  -,  in  the  positive  z  direction  by  a 
distance  Ap.  The  resulting  net  charge  density  in  Figure  4  is  approximated  by 
6-functions. 


p  -  e<n>Ap[6(z)  -  6(z  -  «) 1 
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CHARGE  DENSITY 


DISTANCE  (Z) 


(a)  IRRADIANCE  DISTRIBUTION 

0000 
0000 
0  0  ©  © 

0000 

DISTANCE  (Z)  Apv 

(b)  RESULTING  CHARGE  ARRANGEMENT 


I 


1 

DISTANCE  (Z) 


(c)  RESULTING  CHARGE  DENSITY.^ 

Figure  4.  Model  for  the  photoelectric  effort  showing  the 
idealized  charge  distribution  resulting  from 
net  electron  emission  in  the  positive  z 
direction. 


where  the  electron  concentration  <n>  develops  In  a  characteristic  time  Tr. 
Integrating  Poisson's  equation,  edE/dz  =  4rp  gives, 

E  -  E  for  0  <  z  <  « 

psra 

■  0  for  all  other  values  of  z 

where  EpSm  Is  given  by  Equation  9.9.  Since  Ep  develops  in  time 
^  and  Epsm  develops  in  time  Tr,  It  is  reasonable  that  Epsm 
-  (Tr/Tdl)Ep  *■  (Ydl/Yr>E»  as  ln  Equation  9.9. 
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SECTION  10 


MODEL  OF  DIELECTRIC  RELAXATION 

The  following  model  illustrates  the  significance  of  the  dielectric 
relaxation  frequency: 


=  1/r^  ”  4iro/e  ■  4ne<n>/e  , 


(10.1) 


which  is  the  time  required  to  change  or  discharge  a  sample.  Consider  a 
dielectric  sample  in  an  externally  applied  electric  field  Eapp>  with  open 
circuit  conditions,  as  illustrated  in  Figure  5.  At  time  t  *  0,  the  sample  is 
uncharged.  The  electric  field  Eapp  causes  a  current 


J 


0 


ey<n>E 

app 


at  time  t  **  0.  The  divergence  of  the  current  at  the  end-faces  of  the  sample 
causes  the  end-faces  to  charge  up,  according  to  the  continuity  equation 

3p/3t  =  -3J/3dz  (10.2) 

where  p  is  the  charge  density.  For 

J(z,t)  =•  -ey<n>E(t)r(z)  (10.3) 


where 


r(z)  *  1  for  0  <  z  <  I 

■  0  for  all  other  values  of  z 


(10.4) 


Equation  10.2  gives 


3p/3t  =  -eii<n>E(t)A(z) 


(10.5) 


Figure  5.  Crystal  with  an  applied  field  with 

open-circuit  conditions  at  the  faces 
of  the  crystal  showing  the  current 
density  J,  electric  field  E,  and 
charged  faces  of  the  crystal. 


► 
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« 


1 


where 

A(z)  -  6(z)  -  6(z  -  *)  .  (10.6) 

C 

Here  6(z)  Is  the  Dirac  6-function,  and  dr(z)/dz  =  A(z). 

From  Poisson’s  equation, 

^  P  =  E(4ir)'13E/3t  -  e(d4if)_13[E(t)r(z)]/3z  (10.7) 

*•  ed(4ir)  *E(t)A(z) 


Taking  the  time  derivative  of  p  with  respect  to  time  in  Equation  10.7  and 
substituting  the  result  into  Equation  10.5  gives 

3E(t)/3t  =  -(4ir  ep<n>/e)E(t)  . 


For  E(t)  =  Eapp  at  t  -  0,  the  solution  to  Equation  10.8  is 


(10.8) 


E(t)  -  Ee 
app 


(10.9) 


ft 


M 


where  y^i  is  the  dielectric  relaxation  frequency  defined  in  Equation  10.1. 

The  model  illustrates  that  the  time  constant  for  the  development  of  the 
final  charge  density  (p  *  pgA(z)),  which  corresponds  to  the  charged  end  faces 
of  the  sample)  is  the  dielectric-relaxation  time.  The  electric  field  and  current 
decay  to  zero  as  the  end  faces  of  the  sample  charge  up. 

Using  the  same  model,  but  with  F-app  =  0  and  p  =  pgA(z)  at  time  t  *  0, 
it  is  easy  to  show  that  an  initial  charge  decays  with  time-constant  I/yj^: 


E 


E 


0 


t 


(10.10) 


The  current  and  electric  field  decay  to  zero  as  the  end-face  charge  decays  to 
zero. 
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SECTION  11 


SUMMARY  OF  RESULTS 


A  list  of  symbols  and  typical  values  of  parameters  is  given  In  the  Appendix. 
Key  results  for  the  grating  electric  field  Eg,  which  is  the  z-dependent  part  of 
E  *  ®app  +  ®sc  (which  excludes  the  electric  field  of  the  laser 
irradiance),  are  listed  below. 

The  corresponding  value  of  the  grating  index  is  given  by  Equation  1.1  and 
repeated  here 


n 


g 


1 

2 


n  rE 


(1.1) 


and  the  corresponding  grating  efficiency  Is  given  by  an  expression  such  as 
Equation  1.5  and  repeated  here,  i.e., 


n 


e~ad/cos 


0 


sin 


irn  n  d 
-X.-g  -  . 

2X  cos  9 


(1.5) 


The  notation  for  the  v«r'ous  cases  of  E  below  Is  as  follows: 


Eapp:  applied  field  dominated 


D:  diffusion  dominated 


p:  photovoltaic-effect  dominated 


t:  short-time  solution  (rr  t  t^) 


ss:  steady-state 


"lng- 


ftshf 

ngbl: 


kgA  »  1 
kgA  «  1 

negligible  depletion  and  saturation. 


Combining  the  results  of  Sections  3  through  7  gives  the  following  key 


results: 


From  Equations  3.18,  4.9,  5.7  and  7.7,  for  kgA  «  1 


E  2  yj  te  di  [(E  +  E  )  M  cos  k  z  +  E_,  M  sin  k  z] 

di  LV  app  p'  g  2  d  g 


c/E  K  M  cos  k  z  -  Jrt/eu<n>  +  E  ,  M  sin  k  z' 

(— E - 8 - Q - A _ g_ 

\  _  1  +  M  cos  k  z 

g 

c  (i  -  e  dl  )  .for  Ashort,  negligible 


1  +  M  cos  k  z 
g 


(11.1) 


where 


JQ/ep<n>  “  0 


E 

P 

E  -  E 
P  app 


for  open  circuit 

for  short  circuit 

for  E  with  ohmic  contents 
app 


(11.2a) 

(11.2b) 

(11.2c) 


From  Equations  3. 28,  4.11,  5.9,  and  7.20,  for  kgA  »  1, 


-Y..t 

E  ■  Y,.te  (E.  +  E  )  M  sin  k  z 
di  A  p'  g 


li  +  J±M  .  ,  ,  r  M  74  .  _ 

+  E  - :  “  sin  k  z  + - . - r  sin  2k  z 

.  q  1+irh2  8  1  +  \  r  M2  8 


(11.3) 


x  \1  -  e 


'"din2 
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The  results  for  the  phase  shifts  between  the  grating  electric  field  and  the 
grating  lrradance  <I>  H  cos  kgZ  are  summarized  for  the  various  combination  of 
long-  and  short-  electron-diffusion  distances  and  drift  (in  applied  or 
photovoltaic  field)  and  diffusion-limited  transport  in  Table  2. 


Table  2.  Spatial  Phase  Shift  Between  the  Grating 
(i.e«,  Photof ractive  Electrtic  Field  E)  and  Irradiance  I. 


Drift 


r 


kgA  «  1 


kgA  »  1 


Diffusion 


Drift 

Diffusion 

ir/2 

*/2 

Phase 


-ir/2 


ir/2 
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APPENDIX  B 


LIST  OF  SYMBOLS,  EQUATIONS,  TYPICAL  VALUES  OF  PARAMETERS, 
UNITS,  AND  CONVERSION  FACTORS 


The  values  of  parameters  listed  below  are  only  estimates.  They  are  for 
lithium  niobate.  Values  of  most  of  the  parameters  are  unknown.  The  concen¬ 
trations  C£  and  c+  of  full-  and  empty-optical-donor  levels  are  explicitly 

displayed  in  the  equation.  Two  typical  values  of  interest  are  c  =  0.2  percent 
3 

Fe  =  2  x  10  ppm  of  Fo  for  very  heavily  doped  lithium  niobate  and  <•  =  100  ppm 
for  moderate  doping. 

c  (ppm)  *  concentration  of  N  in  atomic  parts  per  million 
C  6 

c+(ppm)  =  concentration  of  N+  in  atomic  parts  per  million 
cg  or  c+  *  100  ppm,  for  Case  I 

3 

■  2  x  10  ppp  *  0.2  percent,  for  Case  II 

D  *  electron  diffusion  coefficient 

-  pkgT/e 

2 

■  0.46  cm  /sec 


B-l 


■  diffusion  electric  field 


-  Dkg/u  -  kjTkg/« 

■  24.2  statvolts/ca  •  7.7S  kV/ca 

-  photovoltaic  aloe trie  field 


E 

q 


c^(^a)  2.2*  V/ca  -  c^(ppa)  7.6  x  10~3  atatvolts/ca 
characteristic  steady-state,  long-transport-dlstance  electric  field 

4"‘'M-K)/Lkg 

c+(ppa)  204  statvolts/ca  -  c^(ppa)  61  kV/ca 
characteristic  electron-drift  electric  field 

V“k« 

c+(ppa)  1.0  kV/ca 
photon  energy 

4.0  x  10"19  J  -  4.0  x  10"12  erg  -  2.5  eV,  for  X  -  500  no 
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( 


I  ■  assumed  sinusoidal  irradiance  distribution  for  writing  a  grating 

■  <I>  (1  +  M  cos  k  z) 
g 

<I>  *  spatial  average  of  I 

2  9 

=  130  W/cm  =  1.3  x  10  erg/scm,  for  10  mW  at  diameter  d  =  0.1  mm 


J 

P 


photovoltaic  current 


=  pN  I 
r  e 


=  2.8  x  105  cm'1  *  (3.6  x  10~6  cm)'1 
M  =  modulation  factor  in  the  irradiance 

2+ 

Ng  =  concentration  of  Fe  ,  filled  optical-donor  levels 


B-3 


N  n  =  value  of  N  at  1=0. 
el)  e 


16  -3 

c(ppm)  x  5  x  10  cm 


18  -3 

5  x  10  cm  for  c  *  100  ppm 


i n20  -3 

10  cm 


for  c  *  0.2  percent  “  2  x  10 J  ppm 


"ei  5  "e  -  Ne0 


N+  -  concentration  of  Fe  ,  empty  optical-donor  levels,  that  is, 
recombination  levels 


*  value  of  N+  at  I  *  0 


c(ppm)  x  5  x  10 cm-3 


n+a  s  n+  -  N+0 


<n>  =  spatial  average  of  the  conduction-electron  density 


(<YI>/Yr)  Ne()  -  oD  <I>  Ne0/h,„Yr 


^ce(ppm)/c+(ppm)]  (1/130  Wcm'2)  3.2  x  lo' 


9  -3 

cm 


n  =  grating  index 
8 


=  1/2  n  rE 


6.1  x  10  for  E  ®  10  kV 


n  =  real  part  of  the  refractive  index 
r 


=  2.3 


p  =  photovoltaic  coupling  constant  (in  the  expression  J  =  pN  I) 

P  e 


=  eA  CT_./'hw 
P  D 


-  9.6  x  10-25  em7/6/e rg1/2 


r  =  appropriate  linear  electro-optic  coefficient  (in  Equation  1.1) 


=  10~8  cm/V 


from  Glass 


1 


s  =  cross  section  for  recombination,  which  is  the  capture  of 
electron  by  an  empty  optical-donor  level  , 


=  10  cm^,  reasonable  value  for  a  deep,  charged  impurity 


t  =  time  of  laser  irradiation  or  time  to  erase 


v  ^  *  thermal  velocity  of  conduction  electrons 


=  10  cm/s  at  room  temperature 
GREEK  SYMBOLS: 

a  m  optical  absorption  coefficient 


40  cm  *  for  0.2  percent  Fe  in  the  melt  (assumed  to  give 
0.1  percent  Fe  in  the  crystal)* 
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dielectric  constant 


(2. 3) 2  -  5.29 

dielectric  relaxation  frequency 

4iteu  <n>/c 

4ttep  aD  <1>  Ne0/eh<ovths 

£ce(ppm)/c+(ppm)J  (1/130  Wcm  2)  1.93  x  10^  s  * 

|^ce(ppm)/c+(ppm)J  (1/130  Wcm  2)  (5.17  x  10  2  ms)  * 
spatial  average  of  the  optical-pumping  frequency  Yj 

Op  <I>/h,U) 

325  s  *  ■  1/3.1  msec 

1/t^  =  recombination  of  electrons  and  empty  optical-donor 
levels  for  the  case  of  negligible  saturation,  N+  * 

Vth  8  N+0 

9  -1 

c+(ppm)  5  x  10  sec  *  c+(ppm)/0.2  nsec 


6 


electron-diffusion  transport  distance 


(D/Yr)1/2 

c+(ppm)  9.6  x  10  ^  cm 
electron-drift  transport  distance 
ME/yr 

c+(ppm)  ^  3.5  x  10  ^  cm  for  E  =  10  kV/cm 
photovoltaic  characteristic  diffusion  length 


-9  1 

8  x  10  cm,  from  Glass 

wavelength  of  the  laser  radiation 

500  nm 

electron  mobility 

3  2  2 

5.3  x  10  cm  /s  statvolt  =  17.6  cm  /Vs 

laser  frequency 

(500  nm/A)  3.77  x  1015  sec-1 
charge  density 
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optical  cross  section  for  optical  excitation  of  an  electron 
into  the  conduction  band 


-18  -1  1 
10  cm  ,  estimated  from  Glass'  values  of  a  =  40  cm 

for  Ne0  =0.5  (2  x  10"3)  5  x  1022  cm*3;  a D  =  a/NeQ  = 

8  x  10-19  cm”^ 

half  angle  in  the  material  between  the  direction  of 
propagation  of  the  two  beams  that  write  the  grating. 

3 

Kukhtarev-and-coworker  transport  parameter 


(E 

app 


+  *  >/*„ 
p  q 


0. 16/c+(ppm) 


Kukhtarev-and-coworker  transport  parameter 


The  following  electrostatic  units  (esu)  and  conversions  are  listed  for 
the  convenience  of  the  reader.  Brackets  denote  "the  units  of." 

[e]  =  (erg  cm)1^2 

[eu]  =  cm3/sec 

[E]  =  (erg/cm3) 1/2 

[li]  =  cm"^2/erg^2  sec 

[statvolts]  =  (erg/cm)1^2 

1  erg  =  10- 7  J 

1  statamp  =  3.34  x  10~10  A 


1  statvolt  =  300  V 
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